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We review the status of the theory and phenomenology of heavy-quark symmetry, ex- 
clusive weak decays of B mesons, inclusive decay rates and lifetimes of b hadrons, and 
CP violation in S-meson decays. 

1. Introduction 

The rich phenomenology of wealf decays has always been a source of information 
about the nature of elementary particle interactions. A long time ago, /3- and 
/i-dccay experiments revealed the structure of the effective flavour-changing inter- 
actions at low momentum transfer. Today, weak decays of hadrons containing 
heavy quarks are employed for tests of the Standard Model and measurements of 
its parameters. In particular, they offer the most direct way to determine the weak 
mixing angles, to test the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix, and to explore the physics of CP violation. On the other hand, hadronic 
weak decays also serve as a probe of that part of strong-interaction phenomenology 
which is least understood: the confinement of quarks and gluons inside hadrons. 

The structure of weak interactions in the Standard Model is rather simple. 
Flavour-changing decays are mediated by the coupling of the charged current Jqq 
to the W-boson field: 

£cc = -^JScW^i + h.c., (1) 

where 

-^CC = (^6' ^M' ^r) 7^ ML ^ + (ml, Cl, ^l) 7^ ^CKM ( ^L | (2) 

contains the left-handed lepton and quark fields, and 

Vud Vus Vub 

VCKM = I Vcd Vcs Vcb I (3) 

Vtd Vts Vtb, 

is the CKM matrix. At low energies, the charged- current interaction gives rise to 
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local four-fcrmion couplings of the form 
where 

= ^Y7T = 1-16639(2) GeV-^ (5) 

is the Fermi constant. 




Figure 1: Examples of leptonic {B t Vr), semileptonic {B'^ 

D'^e~i>e), and non-leptonic (_B" D+tt^) decays of B mesons. 

According to the structure of the charged-current interaction, weak decays of 
hadrons can be divided into three classes: leptonic decays, in which the quarks 
of the decaying hadron annihilate each other and only leptons appear in the final 
state; semileptonic decays, in which both leptons and hadrons appear in the final 
state; and non-leptonic decays, in which the final state consists of hadrons only. 
Representative examples of these three types of decays are shown in Fig. |[ The 
simple quark-line graphs shown in this figure are a gross oversimplification, however. 
In the real world, quarks are confined inside hadrons, bound by the exchange of soft 
gluons. The simplicity of the weak interactions is overshadowed by the complexity 
of the strong interactions. A complicated interplay between the weak and strong 
forces characterizes the phenomenology of hadronic weak decays. As an example, a 
more realistic picture of a non-leptonic decay is shown in Fig. |^. 

The complexity of strong-interaction effects increases with the number of quarks 
appearing in the final state. Bound-state effects in leptonic decays can be lumped 
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Figure 2: More realistic representation of a non-leptonic decay. 

into a single parameter (a "decay constant"), while those in semileptonic decays 
are described by invariant form factors, depending on the momentum transfer 
between the hadrons. Approximate symmetries of the strong interactions help to 
constrain the properties of these form factors. For non-leptonic decays, on the other 
hand, we are still far from having a quantitative understanding of strong-interaction 
effects even in the simplest decay modes. 

Over the last decade, a lot of information on heavy-quark decays has been col- 
lected in experiments at e'^e~ storage rings operating at the T(4s) resonance, and 
more recently at high-energy e+e~ and hadron colliders. This has led to a rather 
detailed knowledge of the flavour sector of the Standard Model and many of the 
parameters associated with it. There have been several great discoveries in this 
field, such as B^-B^ mixinJ3'i, b ^ u transitionsi^i, and rare decays induced by 
penguin operatorJl'l. Yet there is much more to come. Upgrades of the existing 
facilities at Cornell and Fermilab, as well as the i3-factories to be operated at SLAC, 
KEK, HERA-B and LHC-B, will provide a wealth of new results within the coming 
years. The experimental progress in heavy-flavour physics has been accompanied by 
a significant progress in theory, which was related to the discovery of heavy-quark 
symmetry, the development of the heavy-quark effective theory, and the establish- 
ment of the heavy-quark expansion for inclusive decay rates. The excitement about 
these developments is caused by the fact that they allow (some) model-independent 
predictions in an area in which "progress" in theory often meant nothing more than 
the construction of a new model, which could be used to estimate some strong- 
interaction hadronic matrix elements. In section , we explain the physical picture 
behind heavy-quark symmetry and discuss the construction, as well as simple ap- 
plications, of the heavy-quark effective theory. Section |^ deals with applications of 
these concepts to exclusive weak decays of B mesons. Applications of the heavy- 
quark expansion to the description of inclusive decay rates and lifetimes of h hadrons 
are the topic of section |[ Section || is devoted to a discussion of CP violation in 
meson decays and the physics of the unitarity triangle. 
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2. Heavy-Quark Symmetry 

This section provides an introduction to the ideas of heavy-quark symmetryl§~lll 
and the heavy-quark effective theoryQ~§fl, which provide the modern theoretical 
framework for the description of the properties and decays of hadrons containing a 
heavy quark. For a more detailed_description of this subject, the reader is referred 
to the review articles in Refs.E3~E3. 



2.1. The Physical Picture 

There are several reasons why the strong interactions of systems containing heavy 
quarks are easier to understand than those of systems containing only light quarks. 
The first is asymptotic freedom, the fact that the effective coupling constant of 
QCD becomes weak in processes wi th la rge momentum transfer, corresponding to 
interactions at short-distance scaleSO. At large distances, on the other hand, 
the coupling becomes strong, leading to non-perturbative phenomena such as the 
confinement of quarks and gluons on a length scale -Rhad 1/Aqcd 1 fm, which 
determines the size of hadrons. Roughly speaking, Aqcd 0.2 GeV is the energy 
scale that separates the regions of large and small coupling constant. When the 
mass of a quark Q is much larger than this scale, mg 3> Aqcd, it is called a heavy 
quark. The quarks of the Standard Model fall naturally into two classes: up, down 
and strange are light quarks, whereas charm, bottom and top are heavy quarks? 
For heavy quarks, the effective coupling constant as{mQ) is small, implying that 
on length scales comparable to the Compton wavelength Aq ~ 1/wq the strong 
interactions are perturbative and much like the electromagnetic interactions. In 
fact, the quarkonium systems {QQ), whose size is of order \q / as(rnQ) <C i?had, 
are very much hydrogen-like. Since the discovery of asymptotic freedom, their 
properties could be predictedH before the observation of charmonium, and later of 
bottomonium states. 

Systems composed of a heavy quark and other light constituents are more com- 
plicated. The size of such systems is determined by i?had, and the typical mo- 
menta exchanged between the heavy and light constituents are of order Aqcd • The 
heavy quark is surrounded by a most complicated, strongly interacting cloud of 
light quarks, antiquarks, and gluons. In this case it is the fact that Aq ^ i?had, 
i.e. that the Compton wavelength of the heavy quark is much smaller than the size 
of the hadron, which leads to simplifications. To resolve the quantum numbers of 
the heavy quark would require a hard probe; the soft gluons exchanged between 
the heavy quark and the light constituents can only resolve distances much larger 
than Aq. Therefore, the light degrees of freedom are blind to the flavour (mass) and 
spin orientation of the heavy quark. They experience only its colour field, which 
extends over large distances because of confinement. In the rest frame of the heavy 
quark, it is in fact only the electric colour field that is important; relativistic ef- 



"Ironically, the top quark is of no relevance to my discussion here, since it is too heavy to form 
hadronic bound states before it decays. 
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fects such as colour magnetism vanish as mg oo. Since the heavy-quark spin 
participates in interactions only through such relativistic effects, it decouples. That 
the heavy-quark mass becomes irrelevant can be seen as follows: as mg — > cxd, the 
heavy quark and the hadron that contains it have the same velocity. In the rest 
frame of the hadron, the heavy quark is at rest, too. The wave function of the light 
constituents follows from a solution of the field equations of QCD subject to the 
boundary condition of a static triplet source of colour at the location of the heavy 
quark. This boundary condition is independent of mg, and so is the solution for 
the configuration of the light constituents. 

It follows that, in the limit mg — > oo, hadronic systems which differ only in the 
flavour or spin quantum numbera^of the heavy quark have the same configuration of 
their light degrees of freedom0~t2l. Although this observation still does not allow us 
to calculate what this configuration is, it provides relations between the properties 
of such particles as the heavy mesons B, D, B* and D* , or the heavy baryons A;, and 
Ac (to the extent that corrections to the infinite quark-mass limit are small in these 
systems) . These relations result from some approximate symmetries of the effective 
strong interactions of heavy quarks at low energies. The configuration of light 
degrees of freedom in a hadron containing a single heavy quark with velocity v does 
not change if this quark is replaced by another heavy quark with different flavour 
or spin, but with the same velocity. Both heavy quarks lead to the same static 
colour held. For Nh heavy-quark flavours, there is thus an SU{2iV;i) spin-flavour 
symmetry group, under which the effective strong interactions are invariant. These 
symmetries are in close correspondence to familiar properties of atoms. The flavour 
symmetry is analogous to the fact that different isotopes have the same chemistry, 
since to good approximation the wave function of the electrons is independent of 
the mass of the nucleus. The electrons only see the total nuclear charge. The spin 
symmetry is analogous to the fact that the hyperfine levels in atoms are nearly 
degenerate. The nuclear spin decouples in the limit me/rriN — > 0. 

Heavy-quark symmetry is an approximate symmetry, and corrections arise since 
the quark masses are not infinite. In many respects, it is complementary to chiral 
symmetry, which arises in the opposite limit of small quark masses. There is an 
important distinction, however. Whereas chiral symmetry is a symmetry of the 
QCD Lagrangian in the limit of vanishing quark masses, heavy-quark symmetry 
is not a symmetry of the Lagrangian (not even an approximate one), but rather 
a symmetry of an effective theory, which is a good approximation of QCD in a 
certain kinematic region. It is realized only in systems in which a heavy quark 
interacts predominantly by the exchange of soft gluons. In such systems the heavy 
quark is almost on-shell; its momentum fluctuates around the mass shell by an 
amount of order Aqcd ■ The corresponding fluctuations in the velocity of the heavy 
quark vanish as Aqcd/wq — > 0. The velocity becomes a conserved quantity and 



basis of heavy-quark symmetry are model-independent consequences of QCD in a 
well-defined limit. The symmetry-breaking corrections can, at least in principle, be 




Nevertheless, results derived on the 
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studied in a systematic way. To this end, it is however necessary to recast the QCD 
Lagrangian for a heavy quark, 

CQ = Qiip-mQ)Q, (6) 

into a form suitable for taking the hmit mg — s- oo. 

2.2. Heavy-Quark Effective Theory 

The effects of a very heavy particle often become irrelevant at low energies. It is 
then useful to construct a low-energy effective theory, in which this heavy particle no 
longer appears. Eventually, this effective theory will be easier to deal with than the 
full theory. A familiar example is Fermi's theory of the weak interactions. For the 
description of weak decays of hadrons, the weak interactions can be approximated 
by point-like four-fermion couplings, governed by a dimensionful coupling constant 
Gf- Only at energies much larger than the masses of hadrons can the effects of the 
intermediate vector bosons, W and Z, be resolved. 

The process_Df xemoving the degrees of freedom of a heavy particle involves the 
following stepsa~E3: one first identifies the heavy-particle fields and "integrates 
them out" in the generating functional of the Green functions of the theory. This is 
possible since at low energies the heavy particle does not appear as an external state. 
However, although the action of the full theory is usually a local one, what results 
after this first step is a non-local effective action. The non-locality is related to the 
fact that in the full theory the heavy particle with mass M can appear in virtual 
processes and propagate over a short but finite distance Ax ^ 1/M. Thus, a second 
step is required to obtain a local effective Lagrangian: the non-local effective action 
is rewritten as an infinite series of local terms in an Operator Product Expansion 
Roughly speaking, this corresponds to an expansion in powers of 1/AI. 
It is in this step that the short- and long-distance physics is disentangled. The 
long-distance physics corresponds to interactions at low energies and is the same in 
the full and the effective theory. But short-distance effects arising from quantum 
corrections involving large virtual momenta (of order M) are not reproduced in 
the effective theory, once the heavy particle has been integrated out. In a third 
step, they have to be added in a perturbative way using renormalization-group 
techniques. These short-distance effects lead to a renormalization of the coefhcients 
of the local operators in the effective Lagrangian. An example is the effective 
Lagrangian for non-leptonic weak decays, in which radiative corrections from hard 
gluons with virtual momenta in the range between mw and some renormalization 
scale ^ ~ 1 GeV give^rise to Wilson coefficients, which renormalize the local four- 
fermion interactionsE3~Eil. 

The heavy-quark effective theory (HQET) is constructed to provide a simplified 
description of processes where a h eavy quark interacts with light degrees of freedom 
by the exchange of soft gluonsEJ~0. Clearly, toq is the high-energy scale in this 
case, and Aqcd is the scale of the hadronic physics we are interested in. However, 
a subtlety arises since we want to describe the properties and decays of hadrons 



Heavy- Quark Effective Theory 7 



which contain a heavy quark. Hence, it is not possible to remove the heavy quark 
completely from the effective theory. What is possible, however, is to integrate 
out the "small components" in the full heavy-quark spinor, which describe the 
fluctuations around the mass shell. 

The starting point in the construction of the low-energy effective theory is the 
observation that a very heavy quark bound inside a hadron moves more or less with 
the hadron's velocity v, and is almost on-shell. Its momentum can be written as 

= mQ«^ + k^^ , (7) 

where the components of the so-called residual momentum k are much smaller than 
mg. Note that ?; is a four- velocity, so that = 1. Interactions of the heavy 
quark with light degrees of freedom change the residual momentum by an amount 
of order Afc ~ Aqcd, but the corresponding changes in the heavy-quark velocity 
vanish as Aqco/'TtQ ^0. In this situation, it is appropriate to introduce large- 
and small-component fields hy and Hy by 

hy{x) = e^Q"-" P+ Q{x) , Hy{x) = e™'«''-" P„ Q{x) , (8) 

where P+ and P_ are projection operators defined as 

It follows that 

Q(x) = e-^Q"-" [hyix] + Hy{x)] . (10) 

Because of the projection operators, the new fields satisfy ^ hy = hy and / Hy — 
—Hy. In the rest frame, i.e. for — (1,0,0,0), hy corresponds to the upper two 
components of Q, while Hy corresponds to the lower ones. Whereas hy annihilates 
a heavy quark with velocity v, Hy creates a heavy antiquark with velocity v. 

In terms of the new fields, the QCD Lagrangian (H) for a heavy quark takes the 
form 

Cq = hy iv-Dhy ~ Hy {iv-D + 2mQ) Hy + hy i p±Hy + Hy i p±hy , (11) 

where D'^ = — v ■ D is orthogonal to the heavy-quark velocity: v ■ D± ~ 0. 
In the rest frame, D'^ — (0, -D ) contains the spatial components of the covariant 
derivative. From (p^), it is apparent that hy describes massless degrees of freedom, 
whereas Hy corresponds to fluctuations with twice the heavy-quark mass. These 
are the heavy degrees of freedom that will be eliminated in the construction of the 
effective theory. The fields are mixed by the presence of the third and fourth terms, 
which describe pair creation or annihilation of heavy quarks and antiquarks. As 
shown in the first diagram in Fig. ^, in a virtual process, a heavy quark propagating 
forward in time can turn into an antiquark propagating backward in time, and then 
turn back into a quark. The energy of the intermediate quantum state hhH is 
larger than the energy of the incoming heavy quark by at least 2mQ. Because of 
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this large energy gap, the virtual quantum fluctuation can only propagate over a 
short distance Aa; ~ I/toq. On hadronic scales set by i?had = 1/^QCD, the process 
essentially looks like a local interaction of the form 

h^tp^J—ip^h,, (12) 
2m,Q 

where we have simply replaced the propagator for Hy by l/2mQ. A more correct 
treatment is to integrate out the small-component field Hy , thereby deriving a non- 
local effective action for the large-component field hy, which can then be expanded 
in terms of local operators. Before doing this, let us mention a second type of 
virtual corrections involving pair creation, namely heavy-quark loops. An example 
is shown in the second diagram in Fig. ^. Heavy-quark loops cannot be described 
in terms of the effective fields hy and Hy, since the quark velocities inside a loop 
are not conserved and are in no way related to hadron velocities. However, such 
short-distance processes are proportional to the small coupling constant ^^(toq) 
and can be calculated in perturbation theory. They lead to corrections that are 
added onto the low-energy effective theory in the renormalization procedure. 




Figure 3: Virtual fluctuations involving pair creation of heavy quarks. 
Time flows to the right. 

On a classical level, the heavy degrees of freedom represented by Hy can be 
eliminated using the equation of motion {i p — mg) Q = 0. With ([lO|), this gives 

iphy + iip-2mQ)Hy ^0, (13) 

and multiplying by P± one derives the two equations 

~ ivDhy =ip^Hy , {ivD + 2mQ)Hy = ip^hy. (14) 

The second one can be solved to give 

Hv^T^ ^^^^P±hy, (15) 

2mQ + iv-U 

which shows that the small-component field Hy is indeed of order 1/toq. We can 
now insert this solution into the first equation to obtain the equation of motion 
for hy. It is easy to sec that this equation follows from the non-local effective 
Lagrangian 

Ccs = hyivDhy + hyipj^- — — —ip^hy. (16) 

2mQ + iv-D 

Clearly, the second term corresponds to the first class of virtual processes shown in 
Fig. I 



Heavy- Quark Effective Theory 9 



Because of the phase factor in (|l^, the x dependence of the effective heavy-quark 
field hy is weak. In monientum space, derivatives acting on hy produce powers of 
the residual momentum fc, which is much smaller than mg. Hence, the non-local 
effective Lagrangian (M) allows for a derivative expansion in powers of iD/mg: 



1 



iv ■ D 



2m, 



Q 



ip\hy . 



(17) 



Taking into account that hy contains a P+ projection operator, and using the iden- 
tity 



P+tp^ip^P+=P+ 



P+, 



(18) 



where [iD^,iD'^] ~ ig^G^^ is the gluon field-strength tensor, one finds tha' 



1 
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£off = hyivDhy + hy (i D hy + hy (7^, G^"" hy + 0(1/™^) • (l9) 



2mQ AniQ 
In the limit ttiq —> oo, only the first terms remains: 

'Coc — hy iv- D hy . 



(20) 



This is the effective Lagrangian of the HQET. It gives rise to the Feynman rules 
depicted in Fig. ^ 

■ . . i 1 + ^, 



Figure 4: Feynman rules of the HQET (i,j and a are colour indices). A 
heavy quark is represented by a double line labelled by the velocity v. The 
residual momentum k is defined in (Q). 

Let us take a moment to study the symmetries of this Lagrangian0. Since there 
appear no Dirac matrices, interactions of the heavy quark with gluons leave its spin 
unchanged. Associated with this is an SU(2) symmetry group, under which Coo 
is invariant. The action of this symmetry on the heavy-quark fields becomes most 
transparent in the rest frame, where the generators S*' of SU(2) can be chosen as 

S'^l^^O ; [S\S^]^ie^''^S''. (21) 

Here a* are the Pauli matrices. An infinitesimal SU(2) transformation hy ^ {1 + 
ie- S) hy leaves the Lagrangian invariant: 

SCoo^hy[iv-D,ie- S]hy = 0. (22) 
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Another symmetry of the HQET arises since the mass of the heavy quark does 
not appear in the effective Lagrangian. For heavy quarks moving at the same 
velocity, eq. (|2^) can be extended by writing 

Nh _ 

Coo--Y.~^v^^-^^v- (23) 

1=1 

This is invariant under rotations in flavour space. When combined with the spin 
symmetry, the symmetr y g roup is promoted to SU(2iV/i). This is the heavy-quark 
spin-flavour symmetryll3Ej. Its physical content is that, in the limit mg — > oo, the 
strong interactions of a heavy quark become independent of its mass and spin. 

Consider now the operators appearing at order 1 /vtlq in the effective Lagrangian 
(|l9|). They are easiest to identify in the rest frame. The first operator, 

Okin = 7^ hv {iDi_f K ^ K [iD f K , (24) 

2mQ 2mQ 

is the gauge-covariant extension of the kinetic energy arising from the off-shell resid- 
ual motion of the heavy quark. The second operator is the non-Abelian analogue 
of the Pauli interaction, which describes the colour-magnetic coupling of the heavy- 
quark spin to the gluon field: 

^mae — <7ui/ hy > Jly S'B^ hy . (25) 

4toq mg 

Here S is the spin operator defined in (^l|), and = —^e^^'^G^'^ are the components 
of the colour-magnetic field. The chromo-magnetic interaction is a relativistic effect, 
which scales like l/mg. This is the origin of the heavy-quark spin symmetry. 

2.3. Spectroscopic Implications 

The spin-flavour symmetry leads to many interesting relations between the prop- 
erties of hadrons containing a heavy quark. The most direct consequences concern 
the spectroscopy of such statescd. In the limit toq — > oo, the spin of the heavy 
quark and the total angular momentum j of the light degrees of freedom are sepa- 
rately conserved by the strong interactions. Because of heavy-quark symmetry, the 
dynamics is independent of the spin and mass of the heavy quark. Hadronic states 
can thus be classified by_,the quantum numbers (flavour, spin, parity, etc.) of the 
light degrees of freedomc3. The spin symmetry predicts that, for fixed j ^ 0, there 
is a doublet of degenerate states with total spin J = j ± i. The flavour symmetry 
relates the properties of states with different heavy-quark flavour. 

In general, the mass of a hadron Hq containing a heavy quark Q obeys an 
expansion of the form 

niH = mQ+A + +0{l/ml). (26) 

2mQ ^ 
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The parameter A represents contributions arising from terms in the Lagrangian that 
are independent of the heavy-quark massEJ, whereas the quantity Am^ originates 
from the terms of order l/mg in the effective Lagrangian of the HQET. For the 
ground-state pseudoscalar and vector mesons, one can parametrize the contribu- 
tions from the kinetic energy and the chromo-magnetic interaction in terms of two 
quantities Ai and A2, in such a way thatcS 

Am2 = -Ai + 2 [j(J + 1) - 1] A2. (27) 

The hadronic parameters A, Ai and A2 are independent of mg. They characterize 
the properties of the hght constituents. 

Consider, as a first example, the SU(3) mass splittings for heavy mesons. The 
heavy-quark expansion predicts that 

rriBs - rriBa = A^ - A^ -|- 0{l/mb) , 

rriDs - niD^ = As ~ Ad + 0(1/toc) , (28) 

where we have indicated that the value of the parameter A depends on the flavour 
of the light quark. Thus, to the extent that the charm and bottom quarks can both 
be considered sufficiently heavy, the mass splittings should be similar in the two 
systems. This prediction is confirmed experimentally, sinceEJ 

rriBs -ms^ ^ (90 ± 3) MeV, 

mr,^ -TOD^ = (99± 1) MeV. (29) 

As a second example, consider the spin splittings between the ground-state pseu- 
doscalar (J = 0) and vector {J — \) mesons, which are the members of the spin- 
doublet with i — \- The theory predicts that 

m\, — ra\ = 4A2 + 0{l/mb) , 

ml,, -ml, ^ 4X2 + 0(1 /mc). (30) 
The data are compatible with this: 

m|. - m| ~ 0.49 GeV^ , 

ml, - ml 0.55 GeV^ . (31) 
Assuming that the B system is close to the heavy-quark limit, we obtain the value 

A2~0.12GeV2 (32) 



for one of the hadronic parameters in (^7|). This quantity plays an important role 
in the phenomenology of inclusive decays of heavy hadrons. 

A third example is provided by the mass splittings between the ground-state 
mesons and baryons containing a heavy quark. The HQET predicts that 

"^A6 - rnB ^ Abaryon " Amcson + 0{l/mb) , 

"lA, - m_D = Abaryon " Amcson + 0(1/TOc) . (33) 
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This is again consistent with the experimental results 

mA, - ms = (346 ± 6) MeV , 

mA, - TOD = (416 ± 1) MeV , (34) 

although in this case the data indicate sizeable symmetry-breaking corrections. For 
the mass of the Af, baryon, we have used the value 

TOA, = (5625±6) MeV, (35) 

which, is obtained by averaging the result toa;, = (5639 ± 15) MeV quoted in 
RefE3 with the preliminary value toa^ = (5623 ± 5 ± 4) MeV reported by the 
CDF CollaborationcZl. The dominant correction to the relations ( ^3| ) comes from 
the contribution of the chromo-magnetic interaction to the masses of the heavy 
mesonsf which adds a term 3A2/2toq on the right-hand side. Including this term, 
we obtain the refined prediction that the two quantities 

TOA, -TOB--^ = (312±6) MeV, 
2tob 

TOA, - mo - -—^ = (320 ± 1) MeV (36) 
2m£) 

should be close to each other. This is clearly satisfied by the data. 

The mass formula ( |2^ ) can also be used to derive information on the heavy-quark 
masses from the observed hadron masses. Introducing the "spin-averaged" meson 
masses tob ~ j (to^ + Sm^.) ~ 5.31 GeV and rnn = j {mij + 3to£|.) ~ 1.97 GeV, 
we find that 

rub- nic = {ruB - rjiD) j 1 ~ ^0(1/too)1. (37) 

1^ 2mBmD ^ ) 

Using theoretical estimates for the parameter Ai, which lie in the range@~0 

Ai = -(0.4±0.2) GeV^ (38) 

this relation leads to 

TOfc-TOc = (3.40 ±0.03 ±0.03) GeV, (39) 

where the first error reflects the uncertainty in the value of Ai, and the second 
one takes into account unknown higher-order corrections. As will be discussed in 
section the fact that the difference (to{, — TOc) is determined rather precisely 
becomes important in the analysis of inclusive decays of heavy hadrons. 

For completeness, we note that for the pole masses of the heavy quarks we shall 
adopt the values 

rUb = 4.8 GeV , = 1.4 GeV . (40) 

''Because of the spin symmetry, there is no such contribution to the masses of the Aq baryons. 
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The concept of the pole mass of a heavy quark has been the subject of much dis- 
cussion recently. It has been found that there is an unavoidable ambiguity of order 
Aqcd in any definition of the pole mass extending beyond perturbation theoryiyS. 
Formally, this ambiguity appears as a divergence of the perturbation series, which 
relates the pole mass to a renormalized mass defined at short distances, such as 
the "running mass" rnQimq). As long as we work to a finite order in perturbation 
theory, however, this subtlety can be ignored. The values given above will be used 
in connection with one-loop calculations and thus refer to the one-loop pole masses 
of the heavy quarks. 

3. Exclusive Semileptonic Decays 

Semileptonic decays of B mesons have received a lot of attention in recent years. 
The decay channel B D*iv has the largest branching fraction of all B-meson 
decay modes. From a theoretical point of view, semileptonic decays are simple 
enough to allow for a reliable, quantitative description. The analysis of these decays 
provides much information about the strong forces that bind the quarks and gluons 
into hadrons. Schematically, a semileptonic decay process is shown in Fig. |[ The 
strength of the h ^ c transition vertex is governed by the element Vcb of the CKM 
matrix. The parameters of this matrix are fundamental parameters of the Standard 
Model. A primary goal of the study of semileptonic decays of B mesons is to extract 
with high precision the values of \Vcb\ and \Vub\- In this lecture, we will discuss the 
theoretical basis of such analyses. 



6 y 



D. D^. 



Figure 5: Semileptonic decays of a -B meson. 



3.1. Weak Decay Form Factors 

Heavy-quark symmetry implies relations between the weak decay form factors of 
heavy mesons, which are of particular interest. These relations have been derived 
by Isgur and WiseEJ, generalizing ideas developed by Nussinov and WetzeM, and 
by Voloshin and ShifmarOtj. 

Consider the elastic scattering of a _B meson, B{v) — > B{v'), induced by a vector 
current coupled to the b quark. Before the action of the current, the light degrees 
of freedom in the B meson orbit around the heavy quark, which acts as a static 
source of colour. On average, the b quark and the B meson have the same velocity 
V. The action of the current is to replace instantaneously (at t ~ to) the colour 
source by one moving at a velocity v' , as indicated in Fig. |^. H v — v' , nothing 
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happens; the hght degrees of freedom do not reahze that there was a current acting 
on the heavy quark. If the velocities are different, however, the hght constituents 
suddenly find themselves interacting with a moving colour source. Soft gluons have 
to be exchanged to rearrange them so as to form a B meson moving at velocity 
v'. This rearrangement leads to a form- factor suppression, which reflects the fact 
that as the velocities become more and more different, the probability for an elastic 
transition decreases. The important observation is that, in the limit m?, — + oo, the 
form factor can only depend on the Lorentz boost j — v ■ v' that connects the rest 
frames of the initial- and final-state mesons. Thus, in this limit a dimensionless 
probability function £^{v ■ v') describes the transition. It is called the Isgur-Wise 
functioruj. In the HQET, which provides the appropriate framework for taking the 
limit mil —^ oo, the hadronic matrix element describing the scattering process can 
thus be written as 

— {B{v')\ K^rhv \B{v)) = av ■ v') {v + v'r ■ (41) 

Here, and b^' are the velocity-dependent heavy-quark fields defined in It is 
important that the function ^{v ■ v') does not depend on mt. The factor l/niB on 
the left-hand side compensates for a trivial dependence on the heavy-meson mass 
caused by the relativistic normalization of meson states, which is conventionally 
taken to be 

{B{p')\B{p)) - 2mBv'> i2nfS^p-p') . (42) 



Note that there is no term proportional to {v — v')'^ in (41). This can be seen by 
contracting the matrix element with (v — f')^, which must give zero since = 6„ 
and b^if — byi. 





t<tQ t = tQ t>tQ 

Figure 6: Elastic transition induced by an external heavy-quark current. 

It is more conventional to write the above matrix element in terms of an elastic 
form factor Fci{q'^) depending on the momentum transfer ^ [p — p'Y". 

{B{v')\brb\B{v)) ^ F,M^){p + p'Y , (43) 



where p'^'^ — mBV^'\ Comparing this with (41), we find that 



F,i{q^) = ^{vv') , q^ = ^2ml{v ■ v' - I) . (44) 



Because of current conservation, the elastic form factor is normalized to unity at 
q^ = 0. This condition implies the normalization of the Isgur-Wise function at the 
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kinematic point v ■ v' = 1, i.e. for v ~ v': 

e(l) = l. (45) 

It is in accordanee with the intuitive argument that the probabihty for an elastic 
transition is unity if there is no velocity change. Since for v = v' the daughter 
meson is at rest in the rest frame of the parent meson, the point v ■ v' — I is referred 
to as the zero-recoil limit. 

We can now use the flavour symmetry to replace the b quark in the final-state 
meson by a c quark, thereby turning the B meson into a D meson. Then the 
scattering process turns into a weak decay process. In the infinite mass limit, the 
replacement by/ —^ Cy' is a symmetry transformation, under which the effective 
Lagrangian is invariant. Hence, the matrix element 

^ {D{v')\cy,J^'by\B{v)) =aw'){v + v'r (46) 



is still determined by the same function £^(v-v'). This is interesting, since in general 
the matrix element of a flavour-changing current between two pseudoscalar mesons 
is described by two form factors: 

{D{v')\ C7^6 \B{v)) - U{q^) [p + p'Y - I- {q') (p - p'Y ■ (47) 

Comparing the above two equations, we find that 

2 mB±mD 
f±{q = ^{v-v), 
2y/mBniD 

— TO^ + mjj — 2mB'm]j[v ■ v' ~ 1) . (48) 

Thus, the heavy-quark flavour symmetry relates two a priori independent form 
factors to one and the same function. Moreover, the normalization of the Isgur- 
Wise function a.t v ■ v' — \ now implies a non-trivial normalization of the form 
factors f±{q^) at the point of maximum momentum transfer, (/^j^x = ijnB — rno)^'- 

UiqL.)-^P^- (49) 

The heavy-quark spin symmetry leads to additional relations among weak decay 
form factors. It can be used to relate matrix elements involving vector mesons to 
those involving pseudoscalar mesons. A vector meson with longitudinal polarization 
is related to a pseudoscalar meson by a rotation of the heavy-quark spin. Hence, 
the spin-symmetry transformation cj, — > c^, relates the transition matrix element 
for B ^ D to that for B D* . The result of this transformation isll3: 

^ {D*{v', e)\ Cy,rby \B{v)) = ze'"'"^ el v'^vp ' , 



^mBmo* 

^ {D*{v\e)\cy,-f''robv\B{v)) = e*^ {vv' + 1) - v''' e* ■ v ^{vv') 



^yniBmD* 



(50) 
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where e denotes the polarization vector of the D* meson. Once again, the matrix 
elements are completely described in terms of the Isgur-Wise function. Now this is 
even more remarkable, since in general four form factors, V{q^) for the vector cur- 
rent, and Ai{q'^), i ~ 0,1, 2, for the axial vector current, are required to paiametrize 



these matrix elements. In the heavy-quark limit, they obey the relations 



av-v') = V{q^) = Ao{q') = A,{q') 



2mB'm£)t (v ■ v' — 1) . 



(51) 



Equations §|) and summarize the relations imposed by heavy-quark sym- 
metry on the weak decay form factors describing the semileptonic decay processes 
B ^ D £D and B D*£ v. These relations are model-independent consequences of 
QCD in the limit where mh,mc ^ Aqcd- They play a crucial role in the determi- 
nation of |Vcb|. In terms of the recoil variable w = v-v' , the differential semileptonic 
decay rates in the heavy-quark limit become: 



dw 



487r' 



\Vcb? {tub + mof ml {w^ - 1)^/2 ^^{w) , 



G 



dw 



\Vcb\^ {ttib - ruD')^ y/w^ - 1 {w + i f 



4w rn^B — 2wmBmw 



w 



{rriB ~ mD*y 



e{w). (52) 



These expressions receive symmetry-breaking corrections, since the masses of the 
heavy quarks are not infinitely large. Perturbative corrections of order a"(mQ), 
with Q — b or c, can be calculated order by order in perturbation theory. A more 
difficult task is to control non-perturbative power corrections of order (Aqcd/'tiq)"- 
The HQET provides a systematic framework for analysing these corrections. As 
an example, we have discussed in section the 1/toq corrections to the effective 
Lagrangian. For the more complicated case of weakjiecay form factors, the analysis 
of the l/mq corrections was performed by LukeEj. Later, Falk and the present 
author have also analysed the structure of 1 /itlq corrections for both meson and 
baryon weak decay form factorsc^. We shall not discuss these rather technical issues 
in detail, but only mention the most important result of Luke's analysis. Itconcerns 
the zero- recoil limit, where an Malogue of the AdemoUo-Gatto theoremEJ can be 
proved. This is Luke's theoreniEil, which states that the matrix elements describing 
the leading 1/toq corrections to weak decay a.mplitude s van ish at zero recoil. This 
theorem is valid to all orders in perturbation theoryoHEll. Most importantly, it 
protects the^ -B D*iD decay rate from receiving first-order l/mg corrections at 
zero recoiS. A similar statement is not true for the decay B —> D Iv, however. 
The reason is simple but somewhat subtle. Luke's theorem protects only those 
form factors not multiplied by kinematic factors that vanish for v — v' . By angular 
momentum conservation, the two pseudoscalar mesons in the decay B ^ D £ D must 
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be in a relative p- wave, and hence the amphtude is proportional to the velocity \vd\ 
of the D meson in the i3- meson rest frame. This leads to a factor (w^ — 1) in the 
decay rateL_In such a situation, form factors that are kinematically suppressed can 
contributed. 

3.2. Short-Distance Corrections 

In section , we have discussed the first two steps in the construction of the HQET. 
Integrating out the small components in the heavy-quark fields, a non-local effective 
action was derived, which was then expanded in a series of local operators. The 
effective Lagrangian derived that way correctly reproduces the long-distance physics 
of the full theory. It does not contain the short-distance physics correctly, however. 
The reason is obvious: a heavy quark participates in strong interactions through 
its coupling to gluons. These gluons can be soft or hard, i.e. their virtual momenta 
can be small, of the order of the confinement scale, or large, of the order of the 
heavy-quark mass. But hard gluons can resolve the spin and flavour quantum 
numbers of a heavy quark. Their effects lead to a renormalization of the coefficients 
of the operators in the HQET. A new feature of such short-distance corrections 
is that through the running coupling_constant asimg) they induce a logarithmic 
dependence on the heavy-quark masalj. Fortunately, since q;s(toq) is small, these 
effects can be calculated in perturbation theory. 

Consider, as an example, matrix elements of the vector current V —J^n^Q. 
In QCD this current is partially conserved and needs no renormalizationEl. Its 
matrix elements are free of ultraviolet divergences. Still, these matrix elements 
have a logarithmic dependence on mg from the exchange of hard gluons with virtual 
momenta of the order of the heavy-quark mass. If one goes over to the effective 
theory by taking the limit toq — > cxd, these logarithms diverge. Conseoiiently, 
the vector current in the effective theory does require a renormalizationllS. Its 
matrix elements depend on an arbitrary renormalization scale ^, which separates 
the regions of short- and long-distance physics. If is chosen such that Aqcd ^ 
/i <C rnq, the effective coupling constant in the region between /i and toq is small, 
and perturbation theory can be used to compute the short-distance corrections. 
These corrections have to be added to the matrix elements of the effective theory, 
which contain the long-distance physics below the scale fx. Schematically, then, the 
relation between matrix elements in the full and in the effective theory is 

(V^(mQ))QCD = Co(mQ,M) (14)(/i)>HQET + ^^^^^'^^ (1^i(m)>hqet + ■ • • , (53) 

where we have indicated that matrix elements in the full theory depend on mg, 
whereas matrix elements in the effective theory are mass-independent, but do de- 
pend on the renormalization scale. The Wilson coefficients Ci{mQ, ^) are defined by 
this relation. Order by order in perturbation theory, they can be computed from a 
comparison of the matrix elements in the two theories. Since the effective theory is 
constructed to reproduce correctly the low-energy behaviour of the full theory, this 
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"matching" procedure is independent of any long-distance physics, such as infrared 
singularities, non-perturbative effects, the nature of the external states used in the 
matrix elements, etc. 

The calculation of the coefficient functions in perturbation theory uses the pow- 
erful methods of the renormalization group. It is in principle straightforward, yet 
in practice rather tedious. A comprehensive discussion of most of the existing cal- 
culations of short-distance corrections in the HQET can be found in Ref.E3. 

3.3. Model-Independent Determination of \Vcb\ 

We will now discuss some of the most important applications and tests of the above 
formalism in the context of semileptonic decays of B mesons. A model-independent 
determination of the CKM matrix element \Vcb \ based on heavy-quark symmetry can 
be obtained by measuring the recoil spectrum of D* mesons produced in i? — > D*£ v 
decaysil. In the heavy-quark limit, the differential decay rate for this process has 
been given in (p2|). In order to allow for corrections to that limit, we write 

-i2 



dw 487r3 

Aw rn^ — 2wmB'mD* + rnjyt 



1 



1 {niB - mo^ y 



(54) 



where the hadronic form factor J^{w) coincides with the Isgur-Wise function up 
to symmetry-breaking corrections of order asimg) and Aqcd/"1q. The idea is to 
measure the product |T46|.?-"(w) as a function of w, and to extract |V"cb| from an 
extrapolation of the data to the zero-recoil point w — 1, where the B and the D* 
mesons have a common rest frame. At this kinematic point, heavy-quark symmetry 
helps to calculate the normalization with small and controlled theoretical 

errors. Since the range of w values accessible in this decay is rather small (1 < w < 
1.5), the extrapolation can be done using an expansion around w = 1: 

T{w)=T{l)\l~g^{w-l) + ...]. (55) 

The slope is treated as a fit parameter. 

Measurfiments of the recoil spectrum have been performed first by the ARGUSl^ 
and CLEOEI Collaborations in_j3xperiments operating at the T(4s) resonance, and 
more recently by the ALEPHES and DELPHM Collaborations at LEP. As an ex- 
ample. Fig. 1^ shows the data reported by the CLEO Collaboration. The results 
obtained by the various experimental groups from a linear fit to their data are 
summarized in Table ^. The weighted average of these results is 

|K6|^(1) - (34.6 ±1.7) X 10-3, 

f = 0.82 ±0.09. (56) 
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Figure 7: CLEO data for the product |T46|.F(u'), as extracted from the 
recoil spectrum in i? — > D*£ v decays Eil. The hne shows a hnear fit to the 
data. 

The effect of a positive curvature of the form factor has been investigated by StoneS, 
who finds that the value of \Vcb \ may change by up to +4%. We thus increase 
the above value by (2 ± 2)% and quote the final result as 



\Vcb\T{l) ^ (35.3 ± 1.8) X 10" 



(57) 



In future analyses, the extrapolation to zero recoil should be performed including 
higher-order terms in the expansion ( |55| ) . It can be shown in a model-independent 
way that the sha.pe o f t he form factor is highly constrained by analyticity and uni- 
tarity requirementf£3E3. In particular, the curvature at w = 1 is strongly correlated 



with the slope of the fonn factor. For the value of given in (^6|) , one o 
small positive curvatureES, in agreement with the assumption made in Ref. 



Table 1: Values for |Vcfc| J-'(l) (in units of 10 '^) and g'^ extracted from 
measurements of the recoil spectrum in i? ^ D*£ v decays 





\v.b\m (10-3) 




ARGUS 
CLEO 
ALEPH 
DELPHI 


38.8 ±4.3 ±2.5 
35.1 ± 1.9 ± 2.0 
31.4 ±2.3 ±2.5 
35.0 ± 1.9 ± 2.3 


1.17 ±0.22 ±0.06 
0.84 ±0.12 ±0.08 
0.39 ±0.21 ±0.12 
0.81 ±0.16 ±0.10 



Heavy-quark symmetry implies that the general s 
breaking corrections to the form factor at zero recoil ia 



ructure of the symmetry- 



= r/A 1 + 



A 



QCD 



C2 



^^QCD 



VA (1 + (Jl/mO 



(58) 
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where rjA is a short-distance correction arising from a finite renormahzation of the 
flavour-changing axial current at zero recoil, and (Ji/m^ parametrizes second-order 
(and higher) power corrections. The absence of first-order power corrections at zero 
recoil is a consequence of Lukols theorema. The one-loop expression for tja has 
been known for a long timm'E3EZl: 

^ 1 + ^ ( H^IlUI^ In^-I)^ 0.96 . (59) 

The scale M in the running coupling constant can be fixed by adopting the pre- 
scription of Brodsky, Lepage and Mackenzie (BLM)E3, where it is identified with 
the average virtuality of the gluon in the one-loop diagrams that contribute to rjA ■ 
If as{M) is defined in the MS scheme, the result isE^ M ~ 0.51y/rn^/mb. Several 
estimates of higher-order corrections to rjA have been discussed. A renormalization- 
group resummation of logarithms of the type {as In m;,/mc)", as{as hm-ib/nic)" and 
mc/mb{as In mb/wc)" leads toEl-Enl^Ei tja — 0.985. On the other hand, a resumma- 
tion of "renormalon-chain" contributions of the form /Sp^^a", where /3o = 11— ^nf 
is the first coefficient of the QCD /3-function, gives0 tja — 0.945. Using these par- 
tial resummations to estimate the uncertainty gives tja = 0.965 ± 0.020. Recently, 
Czarnecki has improved this estimate by calculating tja at two-loop ordeiO. His 
result, 

= 0.960 ±0.007, (60) 

is in excellent agreement with the BLM-improved one- loop expression (^9|). Here 
the error is taken to be the size of the two-loop correction. 

The analysis of the power corrections Si/m^ is more difficult, since it cannot rely 
on perturbation theory. Three approaches have been discussed: in the "exclusive 
approach" , all 1 /mg ODerators in the HQET are classified and their matrix elements 
estimated, leading toaL3 (Si/m2 = — (3 ± 2)%; the "inclusive approach" has been 
used to derive the bound (5i/m2 < —3%, and to estimate tha10 (5i/m2 — —(7 ±3)%; 
the "hybrid approach" combines the virtues oLthe former two to obtain a more 
restrictive lower bound on Si/m^- This leads toEa 

= -0.055 ±0.025. (61) 

Combining the above results, adding the theoretical errors linearly to be con- 
servative, gives 

J^(l) = 0.91 ±0.03 (62) 

for the normalization of the hadronic form factor at zero recoil. Thus, the corrections 
to the heavy-quark limit amount to a moderate decrease of the form factor of about 
10%. This can be used to extract from the experimental result ( ^ ) the model- 
independent value 



\Vcb\ = (38.8 ± 2.0exp ± 1.2th) X 10-3 



(63) 
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3.4. Bounds and Predictions for g 

The slope parameter in the expansion of the physical form factor in ( |55| ) differs 
from the slope parameter of the Isgur-Wise function by corrections that violate 
the heavy-quark symmetry. The short-distance corrections have been calculated, 
with the resultlli 

f = g^ + {0.16 ±0.02) + Oil/mg). (64) 

Bjorken has shown that the slope of the Isgur-Wise function is related to the form 
factors of transitions of a ground-state heavy meson into excited states, in which 
the light degrees of freedom carry quantum numbers — or by a sum 
rule which is an expression of quark-hadron duality: in the heavy-quark limit, the 
inclusive sum of the probabilities for decays into hadronic states is equal to the 
probability for the free quark transitioii. If one normalizes the latter probability to 
unity, the sum rule takes the formO^Eil 

+ {w-l) |r{™)(HP + (w + If l-i;2Ml4 +0[{w- If] , 

(65) 

where I, m, n label the radial excitations of states with the same spin-parity quan- 
tum numbers. The terms in the first line on the right-hand side of the sum rule 
correspond to transitions into states containing light constituents with quantum 
numbers = \ ■ The ground state gives a contribution proportional to the Isgur- 
Wise function, and excited states contribute proportionally to analogous functions 
Sy\w). Because at zero recoil these states must be orthogonal to the ground state, 
it follows that ^^'-"(l) = 0, and the corresponding contributions to ( |65|) are actually 
of order {w — if . The contributions in the second line correspond to transitions 
into states with = i"*" or Because of the change in parity, these are p-wave 
transitions. The amplitudes are proportional to the velocity \vj\ — {w^ — if^^ of 
the final state in the rest frame of the initial state, which explains the suppression 
factor (w — 1) in the decay probabilities. The functions Tj{w) are the analogues of 
the Isgur-Wise function for these transitionsE^. Transitions into excited states with 
quantum numbers other than the above proceed via higher partial waves and are 
suppressed by at least a factor {w — if. 

For w = 1, eq. ( |65| ) reduces to the normalization condition for the Isgur-Wise 
function. The Bjorken sum rule is obtained by expanding in powers of {w — 1) 
and keeping terms of first order. Takinginto account the definition of the slope 
parameter, ^'(1) — —g^, one finds thatE2[E3 

^' = i+El-V2^(l)l' + 2El-i;^(l)P>J- (66) 
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Notice that the lower bound is due to the prefactor + of the first term in 
and is of purely kinematic origin. In the analogous sum rule for Aq baryons, this 
factor is absent, and consequently the slope paiameter of the baryon Isgur-Wise 
function is only subject to the trivial constraintBH > 0. 

Voloshin has derived another sum rule involving the form factors for transitions 
into excited states, which is the analogue of Jiie "optical sum rule" for the dipole 
scattering of light in atomic physics. It readaSj 

^^^^^^ + 2E4'2 > (67) 

m n 

where Ej are the excitation energies relative to the mass ttim of the ground-state 
heavy meson. The important point is that this relation can be combined with the 
Bjorken sum rule to obtain an upper bound for the slope parameter g^: 

where Smin denotes the minimum excitation energy. In the quark model, one ex- 
pects'^ that i^miii — "rriM — WQ, and one may use this as an estimate to obtain 
< 0.75. 

The above discussion of the sum rules ignores renormalization effects. Both 
perturbative and non-perturbative corrections to (^6|) and (|68| ) can be incorporated 
using the OPE, where one introduces a momentum scale ^ ~ few x Aqcd large 
enough for asil^-) and power corrections of order (Aqcd/m)" to be small, but other- 
wise as small as possible to suppress contributions from excited statesH. The result 
iso gi^ji„(/i) < g^ < gmax(M)i where the boundary values are shown in Fig. || as a 
function of the scale ^. Assuming that the OPE works down to values fj, ~ 0.8 GeV, 
one obtains rather tight bounds for the slope parameters: 

0.5<e2<0.8, 0.5<f<l.l. (69) 

The allowed region for g^ has been increased in order to account for the unknown 
1/toq corrections in the relation (|^. The experimental result given in ( p6| ) falls 
inside this region. 

These bounds compare well with theoretical predictions for the slope parameters. 
QCD sum rules have been used to calculate the slope of the Isgur-Wiae function. 
The results obtained by various authors are g?-^ 0.84±0.02 (Bagan et al.E3) J].7±0.1 
(NeubertE3), 0.70 ± 0.25 (Blok and Shifmania), and 1.00 ± 0.02 (NarisonS). The 
UKQCD CoUaboratioiLhas presented a lattice calculation of the slope of the form 
factor J-{w), yieldingL2l g^ — 0.9ta3lo.2- We stress that the sum-rule bounds in 
( |69| ) are largely model independent; model calculations in strong disagreement with 
these bounds should be discarded. 



'^Strictly speaking, the lowest excited "state" contributing to the sum rule is D + tt, which has an 
excitation-energy spectrum with a threshold at m,r. However, this spectrum is broad, so that this 
contribution will not invalidate the upper bound for derived here. 
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Figure 8: Bounds for the slope parameter following from the Bjorken 
and Voloshin sum rules. 



3.5. Measurement of B ^ D*lv Form Factors 

If the lepton mass is neglected, the differential decay distributions in B D*£D 
decays can be parametrized by three helicity amplitudes, or equivalently by three 
independent combinations of form factors. It has been suggested thal-a^good choice 
for three such quantities should be inspired by the heavy-quark limitE^O. One thus 
defines a form factor hAi{w), which up to symmetry-breaking corrections coincides 
with the Isgur-Wise function, and two form-factor ratios 



R2{w) = 



{tub + ruD'Y 



1 - 



(70) 



{niB + niD'Y _ 

The relation between w and has been given in (|l|). This definition is such that 
in the heavy-quark limit Ri{w) = R2{w) = 1 independently of w. 

To extract the functions hAi{w), Ri{'w) and R2{w) from experimental data is 
a complicated task. However, HQET-based calculations suggest that the w depen- 
dence of th£_Jorm-factor ratios, which is induced by symmetry-breaking effects, is 
rather mildEll. Moreover, the form factor hAi{w) is expected to have a nearly linear 
shape over the accessible w range. This motivates to introduce three parameters 
g\i, Ri and i?2 by 



where ^(1) 



hAi{w) ~ ^(1) 
Ri{w) ~ Ri, 
R2{w) ~ i?2, 

0.91 ± 0.03 from @. 



1 - g\j_{w - 1) 



(71) 

The CLEO Collaboration has extracted 



these three parameters from an analysis of the angular distributions in _B ^ D* 
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decayg23. The result is: 

= 0.91 ±0.15 ±0.06, 

Ri = 1.18 ±0.30 ±0.12, 

i?2 = 0.71 ±0.22 ± 0.07. (72) 

Using the HQET, one obtains an essentiaUy model-independent prediction for the 
symmetry-breaking corrections to Ri , whereas the corrections to i?2 are somewhat 



model dependent. To good approximationE3 

„ Aa^(mr) A 

Ri~l+ V + 7^ - 1-3 ± 0.1 , 

i?2 =i 1 - K =i 0.8 ± 0.2 , (73) 

2mc 

with K ~ 1 from QCD sum rules0. Here A is the "binding energy" as defined in 
the mass-formula (26). Theoretical calculations E10 as well as phenomenological 
2sS@i 



analysest3E3 suggest that A ~ 0.45-0.65 GeV is the appropriate value to be used in 
one-loop calculations. A quark-model calculation of i?i and i?2 gives results similar 
to the HQET predictionsE3: Ri ~ 1.15 and i?2 — 0.91. The experimental data 
confirm the theoretical prediction that i?i > 1 and i?2 < 1 , although the errors are 
still large. 

There is a model-independent relation between the three parameters determined 
from the analysis of angular distributions and the slope parameter extracted from 
the semileptonic spectrum. It reada3 



6 0(1713 — mo* ) 

The CLEO data give 0.07 ± 0.20 for the difference of the slope parameters on the 
left-hand side, and 0.22 ± 0.18 for the right-hand side. Both values are compatible 
within errors. 

The results of this analysis are very encouraging. Within errors, the experiment 
confirms the HQET predictions, starting to test them at the level of symmetry- 
breaking corrections. 

3.6. Decays to Charmless Final States 

For completeness, we will discuss briefly semileptonic i3-meson decays into charmless 
final states, although heavy-quark symmetry does not help much in the analysis of 
these processes. Recently, the CLEO Collaboration has reported a first signal for the 
exclusive semileptonic decay modes B n £D and B — > piD. The underlying quark 
process for these transitions is b —^ uiD. Thus, these decays provide information 
on the strength of the CKM matrix element Vub- The observed branching fractions 
areS: 

B(B^7r£D) _/ (1-34 ±0.45) X 10-4; ISGW, 
"^'^^"t (1.63 ±0.57) X 10-4; BSW, 
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r (2.28t° j^) X 10-4; ISGW, 
■^^'^^ " \(3.88l}:i^)xl0-4; BSW. 



(75) 



There is a significant model dependence in the simulation of the reconstruction 
efficiencies, for which the models of Isgur et al. (ISGW)il and Bauer et al. (BSW)ii 
have been used. 



Table 2: Values for iKib/Kbl extracted from the CLEO measurement of 
exclusive semileptonic B decays into charmless final states, taking \Vcb\ — 
0.040. An average over the experimental results in (75) is used for all 



except the ISGW and BSW models, where the numbers corresponding to 
these models are used. The first error quoted is experimental, the second 
(when available) is theoretical. 



Method 


Reference 


B -^-kID 


B piv 


QCD sum rules 


Narison 
Ball" 
Khodj. & R 


Lor 

ick]i02 


0.159 ±0.019 ±0.001 
0.105 ±0.013 ±0.011 
0.085 ±0.010 


o.oeet^;;];]^ ± 0.003 

0.094+°:°!^ ±0.016 


Lattice QCD 


UKQCn 
APE 10 


103 
4 


0.103 ±o.oi2j:°:gig 

0.084 ±0.010 ±0.021 




pQCD 


Li & Yu 


105 


0.054 ±0.006 




Quark models 


Bswk 

KC106 
1SGW2 


i 


0.093 ±0.016 
0.088 ±0.011 
0.074 ±0.012 


076+" °" 

U.Um_Q 014 

n 056+" °"^ 

U.UOD_Q 007 

u.u/y_o.oi4 



The theoretical description of these heavy-to-light (h u) decays is more model 
dependent than that for heavy-to-heavy (6 — > c) transitions, because heavy-quark 
symmetry does not help to constrain the relevant hadronic form factors. A vari- 
ety of calculations for such form factors exists, based on QCD sum rules, lattice 
gauge theory, perturbative QCD, or quark models. Table |^ contains a summary 
of values extracted for the ratio |Kif,/Vc6| from a selection of such calculations. 
Some approaches are more consistent than others in that the extracted values are 
compatible for the two decay modes. With few exceptions, the results lie in the 
range 



0.06-0.11 , 



(76) 



which is in good agreement with the measurement of \Vub\ obtained from the end- 
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point region of the lepton speetrum in inclusive semileptonic decaysl 



0.08 ibO.Olcxpi 0.02th. (77) 



cxp 

incl 



Clearly, this is only the first step towards a more reliable determination of \Vub\', 
yet, with the discovery of exclusive b —> u transitions an important milestone has 
been met. Efforts must now concentrate on more reliable methods to determine 
the form factors for heavy-to-light transitions. Some new ideas in t his direction 
have been disc ussed r ecently. They are based on lattice calculationsllHa, analyt- 



transitions 



icity constramtsElJiiij, or variants of the form-factor relations for heavy-to-heavy 



4. Inclusive Decay Rates and Lifetimes 

Inclusive decay rates determine the probability of the decay of a particle into the 
sum of all possible final states with a given set of quantum numbers. An example is 
provided by the inclusive semileptonic decay rate of the B meson, T{B Xc£v), 
where the final state consists of a lepton-neutrino pair accompanied by any number 
of hadrons with total charm-quark number Uc — 1. Here we shall discuss t he t hcq 
retical description of inclusive decays of hadrons containing a heavy quarkE 



From the theoretical point of view, such decays have two advantages: first, bound- 
state effects related to the initial state (such as the "Fermi motion" of the heavy 
quark inside the hadron) can be accounted for in a systematic way using the heavy- 
quark expansion, in much the same way as explained in the previous sections; sec- 
ondly, the fact that the final state consists of a sum over many hadronic channels 
eliminates bound-state effects related to the properties of individual hadrons. This 
second feature is based on a hypothesis known as quark-hadron duality, which is an 
important concept in QCD phenomenology. The assumption of duality is that cross 
sections and decay rates, which are defined in the physical region (i.e. the region 
of time-like momenta), are lajculable in QCD after a "smearing" or "averaging" 



procedure has been appliedtjfil. In semileptonic decays, it is the integration over 
the lepton and neutrino phase space that provides a "smearing" over the invariant 
hadronic mass of the final state (so-called "global" duality). For non-lcptonic de- 
cays, on the other hand, the total hadronic mass is fixed, and it is only the fact 
that one sums over many hadronic states that provides an "averaging" (so-called 
"local" duality) . Clearly, local duality is a stronger assumption than global duality. 
It is important to stress that quark-hadron duality cannot yet be derived from first 
principles, although it is a necessary assumption for many applications of QCD. The 
validit y o f global duality has been tested experimentally using data on hadronic r 



decayaHj. A more formal attempt to address the problem of quark-hadron duality 



can be found in Ref. 



123 



Using the optical theorem, the inclusive decay width of a hadron Hh containing 
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a b quark can be written in the form 

r{Hb^X) = -^2lm{Hi,\T\Ht,), (78) 

where the transition operator T is given by the time-ordered product of two effective 
Lagrangians: 

T = iJd^xT{C,gix),CesiO)}. (79) 

In fact, inserting a complete set of states inside the time-ordered product, we find 
the standard expression 

r(i/fc^X) = ^ ^(27r)4 54(p^_p^)|^^|^^^|^^^|2 (80) 

for the decay rate. For the case of semileptonic and non-leptonic decays, CcS is 
the effective weak Lagrangian given in (H), which in practice is corrected for short- 



distance effect sEH^^IlHS arising from the exchange of gluons with virtuahties between 
mw and mf,. If some quantum numbers of the final states X are specified, the sum 
over intermediate states is restricted appropriately. In the case of the inclusive 
semileptonic decay rate, for instance, the sum would include only those states X 
containing a lepton-neutrino pair. 



e 



bb 




i,G"''b 



Figure 9: Perturbative contributions to the transition operator T (left), 
and the corresponding operators in the OPE (right). The open squares 
represent a four-fermion interaction of the effective Lagrangian CeS, while 
the black circles represent local operators in the OPE. 



In perturbation theory, some contributions to the transition operator are given 
by the two-loop diagrams shown on the left-hand side in Fig. ^. Because of the 
large mass of the b quark, the momenta flowing through the internal lines in these 
diagrams are large. It is thus possible to construct an OPE for the transition 
operator, in which T is represented as a series of local operators containing the 
heavy-quark fields. The operator with the lowest dimension, o? = 3, is bb. It arises 
from integrating over the internal lines in the first diagram shown in the figure. The 
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only gauge-invariant operator with dimension d = 4 is b ip b; however, the equation 
of motion imphes that between physical states this operator can be replaced by 
nibbb. The first operator that is different from bb has dimension d = 5 and contains 
the gluon field, ft is given by bgs'J^^G^'^b. This operator arises from diagrams in 
which a gluon is emitted from one of the internal lines, such as the second diagram 
shown in the figure. 

For dimensional reasons, the matrix elements of such higher-dimensional opera- 
tors are suppressed by inverse powers of the heavy-qu^rk mas s. Thus, any inclusive 
decay rate of a hadron Hb can be written in the for 



luark mas s 

mR-R- 



r{Hb Xf) = -j^^;^ <; 4 {bb)H + 4 -2 + • • • (81) 



1927r-^ 

where the prefactor arises naturally from the loop integrations, are calculable 
coefficient functions (which also contain the relevant CKM matrix elements) de- 
pending on the quantum numbers / of the final state, and {0)h are the (normal- 
ized) forward matrix elements of local operators, for which we use the short-hand 
notation 

{0)H = ^{Hb\0\Hb). (82) 

In the next step, these matrix elements are systematically expanded in powers of 
l/rrib, using the technology of the HQET. Introducing the velocity-dependent fields 
bv of the HQET, where v denotes the velocity of the hadron Hb, one finds^ iil ll3 



{bb), = l- ^'^'''\-fi^'''K oil/ml), 
2mg 

(bgsa^,G^'''b)H - 2fil{Hb) + 0(l/mb) , (83) 
where we have defined the HQET matrix elements 

filiHb) = {Hb{v)\ 6.„ (iDf b, \Hb{v)) , 

f^hiHb) - (^6(^)1 fo^^^M-G'^'fe. \Hb{v)) . (84) 

Here (i-D)^ ~ {iv ■ D)^ — (iD)'^; in the rest frame, this is the square of the operator 
for the spatial momentum of the heavy quark. Inserting these results into (pl|), we 
obtain 



(85) 



It is instructive to understand the appearance of the "kinetic energy" contribution 
fi^, which is the gauge-covariant extension of the square of the 6-quark momen- 
tum inside the heavy hadron. This contribution is the field-theory analogue of the 
Lorentz factor (1 — vj^)^^^ ~ 1 — p^/2ml, in accordance with the fact that the 
lifetime, r = l/F, for a moving particle increases due to time dilation. 
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The main result of the heavy-quark expansion for inclusive decay rates is that 
the free quark decay (i.e. the parton model) provides the first term in a systematic 
l/mb expansion, i.e. 

Tim -> Xf) =. 4{l + 0{l/ml)] . (86) 

For dimensional reasons, the free-quark decay rate is proportional to the fifth power 
of the 6-quark mass. The non-perturbative corrections to this picture, which arise 
from bound-state effects inside the hadron Hb, are suppressed by (at least) two 
powers of the heavy-quark mass, i.e. they are of relative order (AqcD/^^fc)^- Note 
that the absence of first-order power corrections is a simple consequence of the 
equation of motion, as there is no independent gauge-invariant operator of dimension 
d = 4 that could appear in the OPE. The fact that bound-state effects in inclusive 
decays are strongly suppressed explains a posteriori the success of the parton model 
in describing such processes. 

The hadronic matrix elements appearing in the heavy-quark expansion ( |85| ) can 
be determined to some extept from the known masses of heavy hadron states. For 
the B meson, one finds thatEl 

nl{B) = -Ai = (0.4 ±0.2) GeV^, 

^1{B) = 3A2 = ^ {m%, - ml) ~ 0.36 GeV^ , (87) 



where Ai and A2 are the parameters appearing in the mass formula (|2^). For the 
ground-state baryon A5, in which the light constituents have total spin zero, it 
follows that 

MG(Ab) = 0, (88) 
while the matrix element //^(A;,) obeys the relation 

(mA,-mAj-(mB-mz3) = \ fil{B) - fil{Ab)\ ( - ) + 0(l/m|) , (89) 



where Wib and rfi]j denote the spin-averaged masses introduced in connection with 
(^7[). With the value of toa^ given in (^), this leads to 

filiB) - ^2 (^^) ^ (0.01 ± 0.03) GeV^ . (90) 

What remains to be calculated, then, is the coefficient functions for a given 
inclusive decay channel. We shall now discuss the most important applications of 
this general formalism. 

4.1. Determination of \Vcb\ from Inclusive Semileptonic Decays 



The extraction of |V"c&| fr om t 
is based on the expressionlij^ 



r{B XciD) 



10 i nclusive semileptonic decay rate of the B meson 




Glmf. „ I . Ai +3A2 ^ fffll^^ +^£iMla(— 
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6A2 m 



l--f +... , (91) 



where mf, and are the poles mass of the b and c quarks (defined to a given order 
in perturbation theory), and f{x) and g{x) are phase-space functions: 

f{x) = 1- 8a;2 + Sa;*^ - ~ 12x^ In , (92) 



and g{x) is given elsewheret£j. The theoretical uncertainties in this determination 
of \Vcb\ are quite different from those entering the analysis of exclusive decays. In 
particular, in inclusive decays there appear the quark masses rather than the meson 
masses. Moreover, the theoretical description relies on the assumption of global 
quark-hadron duality, which is not necessary for exclusive decays. One should 
distinguish three sources of theoretical uncertainties: 



Hadronic parameters: The non-pcrturbative corrections are very small; with 
-Ai = (0.4 ± 0.2) GeV^ and A2 = 0.12 GeV^, one finds a reduction of the free- 
quark decay rate by —(4.2 ± 0.5)%. The uncertainty in this number is below 1% 
and thus completely negligible. 

Quark- mass dependence: The fact that F ~ ml /{mc/mi,) suggests a strong 
dependence of the decay rate on the value of the 6-quark mass. However, this 
dependence becomes milder if one chooses mf, and Am = — mc as independent 
variables. This is apparent from Fig. |l^, which shows that F ^ mf^ Am^-^. This 
choice of variables is also preferred from a conceptual point of view, since it leads 
to essentially uncorrelated theoretical uncertainties: whereas m^ = — A -|- . . . is 
mainly determined by the A parameter of the HQET, the mass difference Am obeys 
the expansion shown in (|3^), i.e. it is sensitive to the kinetic-energy parameter Ai. 
Theoretical uncertainties of 60 MeV on Am and 200 MeV on mt are reasonable; 
values much smaller than this are probably too optimistic. This leads to 

^ Jf0.10^^V+f0.05^V - • (93) 
Vry„.asscs VV 200 MeV; V 60MeV; ^ ' 



Perturbative corrections: The perturbative corrections are the most su 
part of the analysis. The semileptonic rate is known exactly to ord er as onlyE23, 
although a partial calculation of the coefficient of order existst23. The result is 



2 



= l_1.67^^_(1.68/3o + ...)^^ +.... (94) 

itrec TT V / 

The one-loop correction is moderate; it amounts to about —11%. Of the two- loop 
coefficient, only the part proportional to the /3-function coefficient (3o is known. For 
= 3 light quark flavours, this term gives 1.68/3o — 15.1, corresponding to a rather 
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large correction of about —6%. One may take this as an estimate of the perturbative 
uncertainty. The dependence of the result on the choice of the renoimalization scale 
and scheme has been investigated, too, and found to be of orderlH2l 6%. Yet, the 
actual perturbative uncertainty may be larger than that. A subset of higher-order 
corrections, the so-called "renormalon-chain" contributions of the form /3q can 
be summed to all orders in perturbation theory, leading tclH3 T/Tti-cc — 0.77±0.05, 
which is equivalent to choosing the rather low scale M ~ 1 GeV in This 
estimate is 12% lower than the one-loop result. These considerations show that 
there are substantial perturbative uncertainties in the calculation of the semileptonic 
decay rate. They could only be reduced with a complete two-loop calculation, which 
is however quite a formidable task. At present, we consider ((5r/r)pcrt — 10% a 
reasonable estimate. 

Adding, as previously, the theoretical errors linearly, and taking the square root, 
leads to 

- 10% (95) 

I ych\ 

for the theoretical uncertainty in the determination of jVcfcl from inclusive decays, 
keeping in mind that this method relies in addition on the assumption of global 



quark-hadron duality. Taking the result of Ball et al.ll23 for the central value, we 
quote 



|K.| = (0.0400 ± 0.0040) . (96) 

With the new world averages for the semileptonic branching ratio, Bsl = (10.90 ± 
0.46)% (see below), and for the average B-meson lifetimecEl, tb = (1.60 ± 0.03) ps, 
we obtain 

\Vcb\ = (40.0 ± 0.9exp ± 4.0th) X 10-3 . (97) 



This is in excellent agreement with the value in (63), which has been extracted 
from the analysis of the exclusive decay B — > D*£D. This agreement is gratifying 
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given the differences of the methods used, and it provides an indirect test of global 
quark-hadron duality. Combining the two measurements gives the final result 

|Kb| = 0.039 ± 0.002. (98) 

After Vud and Vus, this is now the third-best known entry in the CKM matrix. 

4.2. Semileptonic Branching Ratio and Charm Counting 

The semileptonic branching ratio of B mesons is defined as 

where Fnl and Frarc are the inclusive rates for non-leptonic and rare decays, respec- 
tively. The main difficulty in calculating _Bsl is not in the semileptonic width, but 
in the non-leptonic one. As mentioned above, the calculation of non-leptonic decays 
in the heavy-quark expansion relies on the strong assumption of local quark-hadron 
duality. 

Measurements of the semileptonic branching ratio have been performed by vari- 
ous experimental groups, using both model-dependent and model- independent anal- 
yses. The status of the results is controversial, as there is a discrepancy between 
low-energy measurements performed at the T(4s) resonance and high-energy mea.- 
surements performed at the resonance. The average value at low energies ial^j 
BsL = (10JlZ.±0.30)%. High-energy measurements performed at LEP, on the other 



hand, givc£2il Bg"^ = (11.11 ± 0.23)%. The superscript (6) indicates that this value 



refers not to the B meson, but to a mixture of b hadrons (approximately 40% 
B^ , 40% B'^, 12% Bs, and 8% A^). Assuming that the corresponding semilep- 
tonic width is close to that of the B mesonf we can correct for this and find 
BsL = (r(S)/r(6)) B^j>^ = (11.30 ± 0.26)%, where t(6) = (1.57 ±0.03) ps is the av- 
erage lifetime corresponding to the above mixture of b hadronscj. The discrepancy 
between the low-energy and high-energy measurements of the semileptonic branch- 
ing ratio is therefore larger than 3 standard deviations. If wc take the average and 
inflate the error to account for this disturbing fact, we obtain 

BsL = (10.90 ±0.46)%. (100) 

In understanding this result, an important aspect is charm counting, i.e. the mea- 
surement of the average number Uc of charm hadrons produced per B decay. Re- 
cently, two new (preliminary) measurements of th is miap tity have been performed. 
The CLEO Collaboration has presented the valueSllsi = 1.16 ± 0.05, and the 
ALEPH Collaboration has reported the resultli23 ric = 1.20 ± 0.08. The average is 



ric = 1.17 ±0.04. 

''Theoretically, this is expected to be a very good approximation. 



(101) 
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In the parton model, one findsE£D Bsl — 13% and ric — 1.15. Whereas Uc is 
in agreement with experiment, the semileptonic branching ratio is predicted to be 
too large. With the establishment of the l/mg expansion the non-pertm'bative 
corrections to the parton model could be computed, and they turned out to be too 
small to improve the prediction. This led Bigi et al. to conclude that values Bsl < 
12.5% cannot be accommodated by theory, ihus giving rise to a puzzle referred to as 
the "baffling semileptonic branching ratio"E2a. The situation has changed recently, 
however, when it was sh own that higher-order perturbative corrections lower the 
value of BgL significantlylnl. The exact order-as corrections to the non-lcptonic 
width have been computed for rric ^ 0, and an analysis of the renormalization 
scale and scheme dependence has been performed. In particular, it turns out that 
radiative corrections increase the partial width T{B — > Xccs) by a large amount. 
This has two effects: it lowers the semileptonic branching ratio, but at the price of 
a higher value of ric- 



The original analysis of Bagan et al. has recently been corrected in an erratum 



139 



Here we shall present the results of an independent nume rica l analysis using the 



same theoretical input (for a detailed discussion, see Ref.E12l). The semileptonic 
branching ratio and Uc depend on the quark-mass ratio rric/rnb and on the ratio 
fi/mb, where fi is the scale used to renormalize the coupling constant as(/i) and 
the Wilson coefhcients appearing in the non-leptonic decay rate. The freedom in 
choosing the scale /i reflects our ignorance of higher-order corrections, which are 
neglected when the perturbative expansion is truncated at order as ■ Below we shall 
consider several choices for the renormalization scale. We allow the pole masses of 
the heavy quarks to vary in the range [see (p9[)] 



TOfc = (4.8±0.2) GeV, - = (3.40 ± 0.06) GeV , (102) 

corresponding to 0.25 < mc/mh < 0.33. Non-perturbative effects appearing at 
order l/m^ in the heavy-quark expansion are described by the single parameter 
A2 — 0.12 GeV'^ defined in (p7|); the dependence on the parameter Ai is the same 
for all inclusive decay rates and cancels out in i?sL and Uc- For the two choices 
fi = mi, and fi — mf,/2, we obtain 



Bsl - 



12.0 ±1.0%; ^l^mb, 
10.9 ±0.9%; 

1.21 T 0.06; fi = mh, 

1.22 T 0.06; ^'^^'^^ 



The uncertainties in the two quantities, which result from the variation of rric/mh in 
the range given above, are anticorrelated. Notice that the semileptonic branching 
ratio has a stronger scale dependence than ric- This is illustrated in Fig. ^ which 
shows the two quantities as a function of /i. By choosing a low renormalization scale, 
values Bsl < 12% can easily be accommodated. The experimental data prefer a 
scale 11/ mi, ^ 0.5, which is indeed not unnatural. Using the BLM scale setting 
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0.29 
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Figure 11: Scale dependence of the theoretical predictions for the semileptonic branching 
ratio and Uc. The bands show the average experimental values given in (|10(]|) and (|101[). 



methodEiLuke et al. have estimated that fi > 0.32to;, is an appropriate scale in 
this casell23. 

The combined theoretical predictions for the semileptonic branching ratio and 
charm counting are shown in Fig. 12, They are compared with the experimental 



results obtained from low- and high-energy measurements. It was argued that the 
combination of a low semileptonic branching ratio a nd a low value of ric would 
constitute a potential problem for the Standard Modelifj. However, with the new 
experimental and theoretical numbers, only for the low-energy measurements a 
small discrepancy remains between theory and experiment. 

Previous attempts to resolve the "problem of the semileptonic branching ratio" 
have focused on four possibilities: 



• It has been argued that the experimental value of Uc may depend on model 
assumptions a b ou t the production of charm hadrons, which are sometimes 



questionabletilrE13. 



• It has been pointed out that the assu rnp tion of local quark-hadron duality 
could fail in non-leptonic B decayslll3~El3. If so, this will most likely happen 
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Figure 12: Combined theoretical predictions for the semileptonic branching ratio and charm 
counting as a function of the quark-mass ratio nic/mb and the renormalization scale /i. The 
data points show the average experimental values for i3sL and ric obtained in low-energy 
(LE) and high-energy (HE) measurements, as discussed in the text. 




in the channel b — > ccs, where the energy release, E — uib — ™x(ccs)j is of 
order 1.5 GeV or less. However, if one assumes that sizeable duality violations 
occur only in this chaimel, it is impossible to improve the agreement between 
theory and experimentlilH. 

Another possibility is that higher-order corrections in the l/mh expansion, 
which were previously thought to be negligible, give a sizeable contribution. 
As will be discussed in more detail below, certain corrections involving the 
participation of a spectator quark are enhanced by phase space, so that they 
lead to effects of relative sizaH] 167r^(AQCD/?Tif,)'^ rather than (AqcD/^T-f,)^. 
They could lower the semileptonic branch ing ratio by up to 1%, depending on 
the size of some hadronic matrix elementallj. Lattice calculations could help 
to confirm or rule out this possibility. 



Finally, there is also the possibility to invoke New PhysicsE_3 tlj. One may. 



for instance, consider extensions of the Standard Model with enhanced flavour- 
changing neutral currents such as 6 — > s 5. The effect of such a contribution 
would be that both Bsl and ric are reduced by a factor (1 + 77 Sg^)^^, where 
V = (Fj-are — F^^g)/FsL- To obtain a sizeable decrease requires values i] ^ 0.5, 
which are large (in the Standard Model, r™g/FsL 0.1), but not excluded 
by current experiments. 

For completeness, we briefly discuss the semileptonic branching ratio for B de- 
cays into a r lepton, which is suppressed by phase space. The ratio of the semilep- 
tonic rates for decays into r leptons and into electrons can be calculated reliably. 
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The result istiS 



B{B ^ Xrur) 
B{B ^Xei?r) 



0.22 ±0.02. 



(104) 



This ratio has been measured at LEP and is found to beE^ 

= 0.234 ± 0.029, (105) 

B{B^XeDr) 
in good agreement with the theoretical prediction. 

4.3. Lifetime Ratios of b Hadrons 

The heavy-quark expansion shows that the lifetimes of all hadrons containing a b 
quark agree up to non-perturbative corrections suppressed by at least two powers 
of 1/to6. In particular, it predicts that 



t{B") 
r{Bd) 
r{B") 



= l + 0{l/ml), 



(1.00 ± 0.01) + 0(l/mg). 



= 1 



f,l{A,)-^l{B) 



2ml 



CG 



f^hiB) 



Oil/ml) 



0.98 + 0(1/771^) 



(106) 



where cq — 1-1, and we have used { pl\ ) and (po|). The uncertainty in the value of 
the ratio t{Bs)/t{B(i) arises from unknown SU(3)-violating effects in the matrix 
elements of B^ mesons. The above theoretical predictions may be compared with 
the average experimental values for the lifetime ratios, which arcO'E£H: 



rjB-) 
t{BO) 

r{Bd) 
r{Ab) 
t{BO) 



1.02 + 0.04, 
1.01 + 0.07, 
0.78 + 0.05. 



(107) 



Whereas the lifetime ratios of the different B mesons are in good agreement with the 
theoretical prediction, the low value of the lifetime of the baryon is surprising. 
To understand the .structure of the lifetime differences requires to go further 



in the 1/mh expansionli£3. Although at first sight it appears that higher-order 



corrections could be safely neglected given the smallness of the l/77i^ corrections, 
this impression is erroneous for two reasons: first, at order l/77i^ in the heavy- 
quark expansion for non-leptonic decay rates there appear four-quark operators, 
whose matrix elements explicitly depend on the flavour of the spectator quark(s) in 
the hadron Hh, and hence are responsible for lifetime differences between hadrons 
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Figure 13: Spectator contributions to the transition operator T (left), 
and the corresponding operators in the OPE (right). Here Ti denotes some 
combination of Dirac and colour matrices. 



with different light constituents; secondly, these spectator effects receive a phase- 
space enhancement factor of IGtt^ with respect to the leading terms in the OPEEIZ. 
This can be seen from Fig. |l^, which shows that the corresponding contributions to 
the transition operator T arise from one- loop rather than two-loop diagrams. The 
presence of this phase-space enhancement factor leads to a peculiar structure of the 
heavy-quark expansion for non-leptonic rates, which may be displayed as follows: 



r = Fo <^ 1 + a;2 
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QCD 

TOfc 



(108) 



Here x„ and y„ are coefficients of order unity. It is conceivable that the terms of 
order IGtt^ (Aqcd/?ti6)'^ could be larger than the ones of order (AqcD/'Tib)^. It is 
thus important to include this type of corrections to all predictions for non-leptonic 
rates. Moreover, there is a challenge to calculate the hadronic matrix elements of 
the corresponding four-quark operators with high accuracy. Lattice calculations 
could help to improve the existing estimates of such matrix elements. 

In total, a set of four four-quark operators is induced by spectator effects. They 
are: 



bj^,{l-j5)taqqj^'il~j5)tab, 
b{l-j5)taqq{i+'-f5)tab, 



rpq 

-'s-p 



(109) 



where q is a light quark, and ta are the generators of colour SU(3). In most previous 
analyses of spectator effects the hadronic matrix e j e rneT its of these operators have 
been estimated making simplifying assumptionsL£3~L£3. For the matrix elements 
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between iJ-meson states the vacuum saturation approximationll^ was assumed, i.e. 
the matrix elements of the four-quark operators have been evaluated by inserting 
the vacuum inside the current products. This leads to 



{B,\T^_^\B,) 



{B,\Ol_^\B,) 
{B,\Tl_p\B,) 



f. 



'■B„ ; 



(110) 



where Jb is the decay constant of the Bq meson, which is defined as 



{Q\qri^b\Bq(v)) ^ifemBy 



(111) 



This approach has been criticized by ChernyakE£3, who estimates that corrections 
to the vacuum saturation approximation can be as large as 50%. 

An unbiased analysis of spectator effects, which avoids assumptions_about had- 



ronic matrix elements, can be performed if instead of (110) one defines 



L4C 



- Bifl 



{Bq\T^_^\Bq) 
{Bg\Tl^p\Bg) 



(112) 



The values of the dimensionless hadronic parameters Bi and are currently not 
known; ultimately, they may be calculated using some field-theoretic approach such 
as lattice gauge theory or QCD sum rules. The vacuum saturation approximation 
corresponds to setting Bi = 1 and = (at some scale /x, where the approximation 
is believed to be valid). For real QCD, however, it is known that 



= 0(1), 



0(l/iVe 



(113) 



where Nc is the number of colours. Below, we shall treat Bi and (renormalized 
at the scale m^) as unknown parameters. Similarly, the relevant hadronic matrix 
elements of the four-quark operators between Ab-baryon states can be parametrized 
by two parameters, B and r, where _B = 1 in the valence-quark approximation, in 
which the colour of the quark fields in the operators is identified with the colour of 
the quarks inside the baryorc 
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4.3.1. Lifetime ratio for B and B'^ 

The lifetimes of the charged and neutral B mesons differ because of two types of 
spectator effects illustrated in Fig. Q They are referred to as Pauli interference 
and W exchangell£J^E£3. In the operator language, these effects are represented by 



the hadronic matrix elements of the local four-quark operators given in (112). In 
fact, the diagrams in Fig. |lj can be obtained from those in Fig. [l3| b y cutting the 
internal lines, which corresponds to taking the imaginary part in (|7^). 
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Figure 14: Pauli interference and W exchange contributions to the hfetimes 
of the B~ and the mesons. The spectator effect in the first diagram arises 
from the interference due to the presence of two identical u quarks in the 
final state. 



The explicit calculation of these spectator effects leads to 



14C 



where 



-0.4Bi + 6.6ei 



-2.2ei + 2.462. 



(114) 



(115) 



Note the factor of IQir^ in (114), which arises from the phase-space enhancement of 
spectator effects. Given that the parton-model result for the total decay width is 



rtot(s) 



3.7x^|y,J^ 



1927r3 

we see that the characteristic scale of spectator contributions is 



47r 



2 fl^B ( '2.71 f. 



nib 



5% 



(116) 



(117) 



Thus, it is natural that the lifetimes of different b hadrons differ by a few per cent. 

The precise value of the lifetime ratio depends crucially on the size of the 
hadronic matrix elements. Taking fs = 200 MeV for the decay constant of the 
B meson (see Ref.E3 and references therein), i.e. absorbing the uncertainty in this 
parameter into the definition of Bi and e^, leads tot 



r{B-) 
t{BO) 



1 + 0.03Bi - 0.71ei + 0.20£2 ■ 



(118) 



The most striking feature of this result is that the coefhcients of the colour-octet 
operators Ty-A and Ts-p are orders of magnitude larger than those of the colour- 
singlet operator Oy-A- As a consequence, the vacuum insertion approximation. 



which was adopted in RefJi£3 to predict that t{B )/T{Bd) is larger than unity by 



an amount of order 5%, cannot be trusted. With of order 1/Nc, it is conceivable 
that the non-factorizable contributions actually dominate the result. Thus, without 
a detailed calculation of the parameters Si no reliable prediction can be obtained. 
Given our present ignorance about the true values of the hadronic matrix elements, 
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we must conclude that even the sign of the sum of the spectator contributions 
cannot be predicted. A lifetime ratio in the range 0.8 < t{B~)/t{B'^) < 1.2 could 
be easily accommodated by theory. 

In view of these considerations, the experimental fact that the lifetime ratio 
turns out to be very close to unity is somewhat of a surprise. It implies a constraint 
on a certain combination of the colour-octet matrix elements, which reads 

£i - 0.3£2 = few %. (119) 



4.3.2. Lifetime ratio for Bg and B^i 



The lifetimes of the two neutral mesons Bg and Bd differ because spectator effects 
depend on the flavour of the light quark, and moreover because the hadronic matrix 
elements in the two cases differ by SU(3) symmetry-breaking corrections. It is 
difhcult to predict the s ign of the net effect, but the magnitude cannot be larger 



than one or two per centEd2Hl^. Hence 



r{Bs) 



1 ± 0{1%) , 



(120) 



which is consistent with the experimental value in (107). Note that t{Bs) denotes 
the average lifetime of the two Bg states. 



4.3.3. Lifetime ratio for Af, and B^ 



Although, as shown in ( 106 ), lifetime differences between heavy mesons and baryons 
start at order l/m1, the main effects are expected to appear at order in the 

heavy-quark expansion. However, here one encounters the problem that the matrix 
elements of four-quark operators are needed between baryon states. Very little is 
known about such matrix elemenla. Bigi et al. have adopted a simple non-relativistic 
quark model and conclude thatll^ 



t{BO) 



0.90-0.95 . 



(121) 



An even smaller lifetime difference has been obtained by Rosnert£2 
An unbiased analysis gives, in the present caselll3: 



t{B°) 



0.98 - 0.17ei + 0.2062 - (0.013 + 0.022S)r , 



(122) 



where B and r are expected to be positive and of order unity. Given the structure 
of this result, it seems difficult to explain the experimental value t{Ai,)/t{B'^) = 
0.78 ± 0.05 without violating the bound (|ll9|)f Essentially the only possibility is to 



■^Another constraint a#i«3s if one does not want to spoil the theoretical prediction for the semilep- 
tonic branching ratio^lj. 
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have r of order 2-4 or so, as there are good theoretical arguments why B cannot 
be much larger than unity. On the other hand, in a constituent quark picture, r is 
the ratio of the wave functions determining the probability to find a light quark at 
the location of the b quark inside the A;, baryon and the B meson, i.e. 



I<''(0)P 



(123) 



and it is hard to see how this ratio could be much different from unity. 

In view of the above discussion, the problem of the short lifetime appears 
as a puzzle, whose explanation may lie beyond the heavy-quark expansion. If the 
current experimental value persists, one may have to question the validity of local 
quark-hadron duality, which is assumed in the theoretical calculation of lifetimes 
and (non-leptonic) inclusive decay rates. 

5. CP Violation 

The violation of CP symmetry is one of the most intriguing aspects of high-energy 
physics. Experimentally, it is one of the least tested properties of the Standard 
Model. To date, there is only a single unambigumas measurement of a CP-violating 



parameter: the measurement of ck in K decaysll£3. The Standard Model description 
of CP violation is very predictive, on the other hand; all CP-violating effects are 
related to the phase S of the CKM matrix. Yet, this description has twomajor 
difhculties: first, CP violation is a necessary prerequisite for baryogenesisli£il, but 
CP violation in the Standard Model is believed to be too small to account for the 
observed baryon asymmetry in the Universe; secondly, there is the so-called strong 
CP problem. The symmetries of the strong interactions allow a term in the QCD 
Lagrangian that violates CP: 

^^TrG^.G^''. (124) 

The problei Ti is th at such a term would contribute to the electric dipole moment of 
the neutror££3ll£3. In general, electric dipole moments of elementary particles are 
sensitive probes of CP-violating effects. Since the only vector that characterizes an 
elementary particle is its spin/ we must have D = dJ. However, J and D have 
different transformation properties under parity and time-reversal transformations. 
Consequently, if either P or T is a good symmetry, we must have d = 0. According 
to the CPT theorem, a violation of time-reversal symmetry implies CP violation. 
This is why measurements of electric dipole moments can be used to constrain CP- 
violating parameters. The current experimental upper bound on the electric dipole 
moment of the neutron 

\dn\ < I -Ix lO^^^e cm (95% CL) (125) 



^In this sense, the term "elementary particle" applies to the neutron, too. 
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implies \6\ < 10 ^, corresponding to an extreme fine-tuning of a parameter in the 
QCD Lagrangian. The above two prohle rns r ,a,1l for extensions of the Standard 



Model, such as the Peccei-Quinn symmetryli£3. The prospects are thus good that 
detailed studies of CP- violating phenomena at future B factories will provide hints 
to physics beyond the Standard Model. 

Besides measurements of the electric dipole moments of the electron and the 
neutron, the most interesting observables for CP violation are the weak decays of 
K and B mesons. In this section, we will present first a general, model-independent 
discussion and classification of CP-violating effects in meson decays. We will dis- 
tinguish three types of CP violation: direct CP violation in weak decays, indirect 
CP violation in the mixing of neutral meson states, and CP violation in the inter- 
ference of mixing and decay. We will then focus on an analysis of these effects in 
the Standard Model. A more detailed discussion of CP violation in and bevaud the 



Standard Model can be found in the comprehensive review articles by Nirll^, and 
by Buchalla et al. 
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5.1. P, C, and CP Transformations 

We start with a discussion of the parity, charge-conjugation and CP transformations 
acting on meson states. The parity transformation is a space-time transformation, 
under which < — > i, x — > —x. It changes the sign of momenta, p — p, leaving spins 
unchanged. For pseudoscalar mesons P and P, the parity transformation implies 
(adopting the common phase conventions) 

P |P(p)) - -|P(-p)> , P |P(p)) = -|P(-P)) . (126) 

Charge conjugation is a transformation that relates particles and antiparticles, 
leaving all space-time coordinates unchanged, i.e. 

C \P{p)) = \P{p)) , C \P{p)) = |P(p)) . (127) 

The combined transformation, CP, acts on the pseudoscalar meson states as 
follows: 

CP |P(p)) = -\Pi-p)) , CP |P(p)) = -\Pi-p)) . (128) 
For neutral mesons, P° and P", one can construct the CP eigenstates 

|P°) = i= (|P") - |P")) , \P^) = (|P°) + |P")) , (129) 

which obey 

CP |P°) - |P°) , CP |P°> - -|P°> . (130) 



5.2. Direct CP Violation in Weak Decays 

Consider two decay processes related to each other by a CP transformation. Let 
P and P be CP-conjugated pseudoscalar meson states, and / and / some CP- 
conjugated final states: 

CP |P) = e^-^^ |P) , CP I/) = e^-^^ I/) . (131) 
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The phases tpp and tpj are arbitrary. The CP-conjugated decay amplitudes, A and 
A, can then be written as 



A = {f\n\P) = Y,A,e''^e^' , 

i 

A = {f\n \P) = e^(^^-^^) e-*-^* 



(132) 



where Ti is the effective Hamiltonian for weak decays, and Ai are real partial ampli- 
tudes. Two types of phases may appear in the decay amplitudes: the weak phases 

are parameters of the Lagrangian that violate CP. They usually appear in the 
electroweak sector of the theory and enter A and A with opposite signs. The strong 
phases Si appear in scattering amplitudes even if the Lagrangian is CP invariant. 
They usually arise from rescattering effects due to the strong interactions and enter 
A and A with the same sign. 

Although the definition of strong and weak phases is to a large degree convention 
dependent, one can show that the ratio 



A 




A 





Aj e'^^ e'-^^ 
A, e'S, e- 



(133) 



is independent of phase conventions and therefore physically meaningful. The con- 
dition 

— 7^ 1 =^ direct CP violation (134) 

implies CP violation, which results from the interference of decay amplitudes leading 
to the same final state. Note that this requires at least two partial amplitudes that 
differ in both their weak and strong phases. 



Experimental observation of direct CP violation: Since mixing is unavoid- 
able in neutral meson decays, it is best to observe direct CP violation in the decays 
of charged mesons. One defines the CP asymmetry: 

_ r(p+^/)-r(p-^/) _ i~|AMp 
^ r(p+ ^ /) + r(p- ^ /) 1 + \A/A\^ ■ ^ ' 

The requirement of at least two partial amplitudes with different phases forces us 
to consider non-leptonic decays, since leptonic and semileptonic decays are usually 
dominated by a single diagram. Non-leptonic decay's, on the other hand, can re- 
ceive so-called "tree" and "penguin" ^contributionsE^. Penguin diagrams contain a 
W^- boson-quark loop and typically involve other weak phases than tree diagrams. 
In order to get large interference effects, one needs partial amplitudes with similar 

"It is a challenge to draw a penguin diagram in such a way that it would actually deserve its 
name. Among the various rumours about the origin of the name "penguin" , the author tends to 
believe the one of a well-known CERN theorist, who had a bet that he could introduce any name 
he wanted into high-energy physics. 
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magnitudelisi. A possibility is to consider decays in which the tree contribution is 
suppressed, with respect to the penguin contribution, by small CKM parameters. 
This compensates for the loop suppression of penguin diagrams. In the Standard 
Model, an example of this type is the decay K^p^ shown in Fig. for which 

the tree diagram is proportional to the small CKM parameters ~ 10"'^, 

whereas the penguin diagram is proportional to (Q!s/127r) ln(m( /m^) iVtfcVj 



tsl 



0.02 X 0.04 ~ 10~'^. Another possibility is to consider tree-forbidden decays, which 
can only proceed through penguin diagrams. In this case, it is the possibility to 
have different quarks in the loop {t, c, u) that leads to the interference. Examples 
are K^K and B^ K^4>^ as well as the radiative decays B^ — > K*^'-f and 

B^ P^7, see Fig. |l^. Unfortunately, the decays — tt^tt" and B^ tt^tt" 
are pure AI = | transitions and are thus governed by a single strong phase 62 
corresponding to a tttt final state with isospin 1 = 2; isospin / = 1 is not allowed 
because of Bose symmetry. It follows that Utt^ = 0. There is no unambiguous 
experimental evidence for direct CP violation yet. 




Figure 15: Tree and penguin diagrams for the decay B^ K^p^. 



1 




Figure 16: Penguin diagrams for some tree-forbidden B decays. 



Hadronic uncertainties: Calculations of direct CP-violating asymmetries have 
large theoretical uncertainties. They are limited by our incapacity to calculate 
hadronic matrix elements of quark operators with high accuracy. Moreover, direct 
CP violation requires non-trivial strong phase shifts, which are notoriously hard to 
calculate. In some cases, however, part of the uncertainty can be eliminated using 



isospin analysis 
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5.3. Indirect CP Violation in the Mixing of Neutral Mesons 

The neutral mesons and P*^ can mix via common decay channels: 

P° ^ X ^ . (136) 
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An arbitrary neutral meson state can thus be written as a superposition of the 
flavour eigenstates, a|P°) + b\P°), which obeys the time-dependent Schrodinger 
equation 



dt 



H 



).(M 



(137) 



where M and T are Hermitian 2x2 matrices, which are called the mass and 
decay matrices, respectively. Since the Hamilton operator, H, is not Hermitian, 
its eigenvectors 

|Pi,2)=p|P°)±?|P°); \p^\ + \q^\ = l (138) 
are not orthogonal, and the eigenvalues 



Mi 



Mi--ri: i=l,2 
2 



(139) 



are complex. This reflects that the states Pi and P2 are resonances, not elementary 
particles. Mj are the masses of these resonances, and Fj are their decay widths. 
The states P, have a diagonal time evolution given by 



One can show that the ratio 



Pirn- 



M* — IF* 

M12 - i F12 



(140) 



(141) 



is independent of phase conventions and therefore physically meaningful. The con- 
dition 



p 



indirect CP violation 



(142) 



implies CP violation, which results from the fact that the flavour eigenstates are 
different from the CP eigenstates. 

Let us collect some useful equations related to the mixing of neutral mesons. 
Define the mass diff'erence Am = 7712 — mi and the width difference AF = F2 — Fi. 
Then the following relations hold: 

(Am)2-i(AF)2 = 4|Mi2|2 - |Fi2^ 

Am- AF = 4Re(Mi2rj2), 



p 



1 Am - f AF 



12 



2 ^ 12 



2 M12 - I F12 

An alternative common notation is to define e such that 



Am - I AF 



(143) 
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q/p in the kaon system: One defines the "short-hved" and "long-lived" neutral 
kaon states Ks = Ki and — K2, which differ significantly in their lifetimes: 
Ts = (8.926 ± 0.012) X 10"" s and tl = (5.17 ± 0.04) x 10"^ s. Experimentally, 



ArriK = rriL-ms^ (3.510 ± 0.018) x lO'^^ GeV , 
ATk = Tl- Ts = -(7.361 ± 0.010) x 10"^^ GeV , 



so that 



If we define 



^ 12 

Ml 



ri2 



M12 



(145) 
(146) 
(147) 



the experimental observation that there is only a small CP violation in the kaon 
system is reflected in the fact that — 0(10^^). From (143), we find to first 

order in this small angle 



9 

P/ K 



l + j( AT^) 

V 2 Am w ' 



1 - 1012 



2AmK - 



l + (,Ar^^2 



(148) 



^ 2AmA- 



so that with (|146|) 



- 1 ~ -2 Re ex =i -(/'12 = 0(10-^) . 



(149) 



K 



q/p in the i3-meson system: Decay channels common to and , which 
are responsible for the difference AF^, are known to have branching fractions of 
order 10~^ or less. Hence, although ATb has not yet been measured directljf 
it follows that |AFb|/Fb < 10"^. The observed mixing rate implieila 

Ams/Fs = 0.74 ± 0.04, on the other hand, so that model independently 



|AFb| < Ams 



(150) 



Thus, there is a negligible lifetime difference between the CP eigenstates, and one 
therefore refers to these states as "light" and "heavy" , B^ = Bi and B^ = B2- It 



follows that IF12I <C IM12I, and to first order in F12/A/12 we obtain from (143) 



Hence 



PJ B 



M 



12 



121 



1 Im — ^— 

2 Afi2 



1 ~ -2Reei3 = 0(10" 



(151) 



(152) 



As in the kaon system, CP violation in B^-B'^ mixing is a small effect. 
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Experimental observation of indirect CP violation in the kaon system: 

One uses the fact that the semileptonic decays of neutral mesons are flavour-tagging, 
i.e. P'^-f^i^i' and P'^-/^£'^D, and defines 

A= {l+DX\n\P°) , A* = {riyXimP") . (153) 

Because of the large lifetime difference between the two neutral kaon states, it is 
possible to prepare a beam of Kl particles and to measure the asymmetry 

r{KL -> £+DX) + T{Kl -^(.-vX)' ^ ' 

Using that \Kl) = p\K^) - q\K^), and hence 

{t+vX\n\KL) =pA, {l-vX\n\KL) = -qA* , (155) 

we find that 



^ _ l-\qlp? _ 2RelK 



2Reeif. (156) 



1 + k/pP l + \lK? 

Experimentally, it was found that 

a^L = (3-27 ±0.12) X 10"^ (157) 
This was the observation of indirect CP violation in the kaon systemE 
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Experimental observation of indirect CP violation in the B-meson sys- 
tem: Since and Bh have almost identical lifetimes, it is not possible to pro- 
duce selectively beams of Bl or Bh particles. With mn.L — rriB ± ^Ato^ and 
Th,l — Tb, eq. (|l40| ) gives for the time evolution of an initially pure state: 



150(0)) ^1^(\Bh) + \Bi 

\B°{t)) ^ — e-™«* e-^rs* J 6-5^"^* \Bh) + e*'^'"^* \Bl) 



2p 

^"^^ e-^^^* <! cos f i ATOst) 15°) + ^ sin f i A77iBi] I^O) 



(158) 



Similarly: 



\B\t)) = e-"^«* e~^^^* I cos (i Amst) \B°) + ^ sin (i Amst) |B0>| . (159) 



Defining the semileptonic asymmetry as 



B ^ r(^"(t) -> l+DX) ~ T{B^{t) -> l-vX) 

T{B^{t)^l+DX)+T{B^{t)^l-vX)' ^ ' 
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and taking into accomit that B^-f^t v and B^-f^l'^v, we obtain 

«fL= I'j^^^j' c.4Rce-B = O(10-^). (161) 

To date, there is no experimental evidence for indirect CP violation in the _B-meson 
system. 

Hadronic uncertainties: The calculation of \q/p\ involves hadronic matrix el- 
ements of local four-quark operators (so-called B parameters). The theoretical 
uncertainty in the calculation of such matrix elements is about 30%. 

5.4. CP Violation in the Interference of Mixing and Decay 

Consider decays of neutral mesons into CP eigenstates: 

^ = (/cp I H |P°) , ^* = (/cP I H |P°) . (162) 
It can be shown that the product 

A = ^ • 4 (163) 

is independent of phase conventions and thus physically meaningful. In other words, 
the convention dependence oi q/p cancels against that of A/ A. The condition 

A 7^ 1 CP violation (164) 

implies CP violation. Note that direct CP violation (|^/^| 7^ 1) and indirect CP 
violation {\q/p\ ^ 1) imply |A| ^ 1, but they are not necessary for the weaker 
condition A 7^ 1. In fact, the case |A| = 1 but ImA 7^ is the theoretically favoured 
situation. In that case A is a pure phase, which can be calculated without hadronic 
uncertainties. 

Many decays of_ncu1jral B mesons are of the kind described above. If one defines 
the CP asymmetr; 



it neu tral 
y&0 



„ r(i3"(t)^/cp)-r(^°(0^fcp) 

r(s"(<)^/cp) + r(BO(t)^/cp) ^ ^ 

and takes into account that \q/p\B — 1, it follows that 

(1 - |Ap) cos(AmB<) - 2ImA sin(AmBt) 



1 + IAp 

'^-i^^ -ImA sin(Amst) . (166) 



The "clean modes" with |A| ~ 1 are those dominated by a single weak phase 0, so 
that 

^ ~ e-^'^ (167) 
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is close to a pure phase. Examples of such decays are discussed in detail below. 
Unfortunately, this method is not useful in kaon decays, since 

Im X{K tttt) ^ 0(10^3) , (168) 

i.e. very small. 

5.5. CP Violation in the Standard Model 

We will now specify the general framework described above and discuss CP violation 
in the context of the Standard Model. Below mass scales of order mw ~ 80 GeV, the 
Standard Model gauge group SUc(3) x SUl(2) x Uy(1) is spontaneously broken to 
SUc(3) X Uom(l), since the scalar Higgs doublet (f) acquires a vacuum expectation 
value. This gives masses to the W and Z bosons, as well as to the quarks and 
leptons. The quark masses arise from the Yukawa couplings to the Higgs doublet, 
which in the unbroken theory are assumed to be of the most general form invariant 
under local gauge transformations. The Yukawa interactions are written in terms 
of the weak eigenstates q' of the quark fields, which have simple transformation 
properties under SUl(2) x Uy(1). After the symmetry breaking, one redefines the 
quark fields so as to obtain the mass terms in the canonical form. This has an 
interesting effect on the form of the flavour-changing charged-current interactions. 
In the weak basis, these interactions have the form 



A^t = {ul cli'^) r\s'^\ + h.c. (169) 



V2 

In terms of the mass eigenstates q, however, this becomes 




Ant = (^L, CL, tL) 7^ ^CKM | SL | + h.C. (170) 

The CKM mixing matrix Vckm is a unitary matrix in flavour space. In the general 
case of n quark generations, Vckm would be an n x n matrix. 

For two generations, the mixing matrix can be parametrized by one angle and 
three phases: 

cos 0c e'" sin 0c e*'' 
-sin^ce*'' cos6'ce^^'3+'^""' 



V = 



e'" \ / cos6lc sin^c \ fl 
e*''M-sin6lc cos^c M e*^''-"' 



(171) 



The phases are not observable, however, as they can be absorbed into a redefinition 
of the phases of the quark fields ml, cl, sl relative to dh- After this redefinition. 



the matrix takes the standard form due to CabibbollZZI: 



cos 0c sm0c , (,72) 
' — sm Oq cos Oc ' 
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In the case of three generations, Vckm can be parametrized by three Euler angles 
and six phases, five of which can be removed by adjusting the relative phases of the 
left-handed quark fields. Hence, three angles 6*^ and one observable phase 5 rer nain 
in the quark mixing matrix, as was first pointed out by Kobayashi and MaskawalilS. 
For completeness, we note that in the general case of n generations, it is e a,sy to show 



that there are ^n(n— 1) angles and \ {n — l)(n — 2) observable phasesliI3. Therefore, 
whereas the Cabibbo matrix is real and has only one parameter, the CKM matrix 
is complex and can be parametrized by four quantities. The imaginary part of the 
mixing matrix is necessary to describe CP violation in the Standard Model. In 
general, CP is violated in flavour-changing decays if there is no degeneracy of any 
two quark masses, and if the quantity Jcp ^ 0, where 

JcP = I Im {V,,VmV*iV^^) I ; i^k, j^l. (173) 

It can be shown that all CP-violating amplitudes in the Standard Model are pro- 
portional to J r'p, a,nd that this quantity is invariant under phase redefinitions of 
the quark fields 



18C 



Ignoring the strong CP problem, i.e. assuming that 6* = in (124), the complex 
phase of the CKM matrix is the only parameter in the Standard Model that violates 
CP symmetry. Hence, the Standard Model is very predictive in describing CP- 
violating effects: all CP-violating observables are in principle calculable in terms 
of only one parameter. In practice, however, strong interaction effects have to be 
controlled before such calculations can be performed. 

Let us mention two of the most convenient parametrizations of the CKM matrix. 



The "standard parametrization"ll£^ recommended by the Particle Data Group is 

Cl2 Cl3 Sl2 Cl3 Sl3 e""' 

Vckm = | -si2 C23 - C12 S23 S13 e"' C12 C23 - S12 S23 S13 e*'' S23 C13 | . (174) 
S12 S23 - C12 C23 si3 e"' -C12 S23 - S12 C23 si3 e"^ C23 ci3 

Here, the short-hand notation Cij = cos Oij and Sij — sin 9ij is used. Some advan- 
tages of this parametrization are the following ones: 

• \Vub\ = Si3 is given by a single angle, which experimentally turns out to be 
very small. 

• Because of this, several other entries are given by single angles to an accuracy 
of better than four digits. They are: Vud — ci2, Vus — 512, Vcb — S23, and 
Vtb — C23. 

• The CP-violating phase S appears together with the small parameter S13, 
making it explicit that CP violation in the Standard Model is a small effect. 
Indeed, one finds 

Jcp = I si3 S23 S12 ss cfg C23 C12 | . (175) 



For many purposes and applications, it is more convenient to use an approximate 
parametrization of the CKM matrix, which makes explicit the strong hierarchy 
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observed experimentally. Setting C13 = 1 (experimentally, it is known that C13 > 
0.99998) and neglecting S13 compared with terms of order miity, we find 



C12 

VcKM — I — S12 C23 

S12 S23 - C12 C23 si3 e 



iS 



S12 
C12 C23 
- C12 S23 



si3 e 
S23 
C23 



(176) 



Now denote A = S12 — 0.22. Experiments indicate that S23 = 0{X^) and S13 = 
O(A^). Hence, it is natural to define S23 — AX'^ and S13 e"'* = A A'^(/9 — ir]), with 
A, p and r] of or der unity. An expansion in powers of A then leads to the Wolfenstein 
parametrizatior 



1 - 



VcKM 



A 



1 - — 



\AX^{1- p~i7]) 



AX^ 



AX^p-iT^)\ 
AX'^ 
1 / 



0{X'). 



(177) 



It nicely exhibits the hierarchy of the mixing matrix: the entries on the diagonal 
are close to unity, Vus and Vcd are of order 20%, Vcb and Vts are of order 4%, and 
Vub and Vtd are of order 1% and thus the smallest entries in the matrix. Some care 
has to be taken when one wants to calculate the quantity Jcp in the Wolfenstein 
parametrization, since the result is of order A^ and thus beyond the accuracy of the 



approximation. However, taking i = u, j = d, k = t, and ^ = 6 in (173), we obtain 
the correct answer 

Jcp ~ A^T^A^ ~ 1.1 X lO^'^A^T?, (178) 

which shows that Jcp is generically of order 10^"* for A ~ 0.22. 

In principle, the entries in the first two rows of the mixing matrix are accessible 
in so-called direct (tree-level) processes, i.e. in weak decays of hadrons containing 
the corresponding quarks. In practice, \Vud\ and \Vus\ are known to an accuracy of 
better than 1%, \Vcb\ is known to 5%, and \Vcd\ and \Vcs\ are known to about 10- 
20%. Hence, the two Wolfenstein parameters A and A are rather well determined 
experimentally: 



X = IK. 



0.2205 ±0.0018. 



A 



y2 



0.80 ±0.04. 



(179) 



On the other hand, \Vub\ has an uncertainty of about 30%, and the same is true 
for I Vtd I, which is obtained from B^-B^ mixing. This implies a rather significant 
uncertainty in the values of the Wolfenstein parameters p and r/. A more precise 
determination of these parameters will be a challenge to experiments and theory 
over the next decade. 



5.6. The Unitarity Triangle 

A simple but beautiful way to visualize the implications of unitarity is provided by 
the so-called unitarity trianglell23, which uses the fact that the unitarity equation 



V, V*^ = (j ^ k) 



(180) 
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can be represented as the equation of a closed trianglciji^the complex plane. There 
are six such triangles, all of which have the same area 



184 



I^aI = \Jcp 



(181) 



Under phase reparametrizations of the quark fields, the triangles change their ori- 
entation in the complex plane, but their shape remains unaffected. 

Most useful from the phenomenological point of view is the triangle relation 



(182) 



since it contains the most poorl y known entries in the CKM matrix. It has been 
widely discussed in the literaturell£l~ll£§. In the standard parametrization, VcdV*f, 
is real, and the unitarity triangle has the form shown in Fig. It is useful to 
rescale the triangle by dividing all sides by y ^^V*^^ The rescaled triangle has the 



coordinates (0,0), (1,0), and (p, 77), wheret^l 

A2 



P= 1- 



^= 1- 



V 



(183) 



are related to the Wolfenstein parameters p and 77 appearing in (177). Unitarity 
amounts to the statement that the triangle is closed, and CP is violated when the 
area of the triangle does not vanish, i.e. when all the angles are different from zero. 



td^tb 





VcdVcb (0.0) (1.0) 

Figure 17: The unitarity triangle (left), and its rescaled form in the p- 
fj plane (right). The angle 7 coincides with the phase S of the standard 
parametrization. 

To determine the shape of the triangle, one can aim for measurements of the two 
sides Rb and Rt, and of the angles a, f3, and 7. So far, experimental information is 
available only on the sides of the triangle. The current value of \Vub\ 



Rb = VW+ 



1 



A2 
2 



Vcl 



0.35 ±0.09. 



in (77) implies 
(184) 



To determine i?^, one needs information on |Vtd|i which can be extracted from B^- 
mixing. In the Standard Model, the mass difference Atob between the two 
neutral meson states is calculable from the box diagrams shown in Fig. |l^. The 
resulting theoretical expression is 



67r2 



^ nB Bsj^ms S{mtlrtiw) iV^dV^^P . 



(185) 
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194 



where rjs — 0.55 ± 0.01 accounts fQi,|tliij QCD correctionsli23, and S{mt/mw) is 
a function of the top quark mass 



The product Bb/b parametrizes the 
hadronic matrix element of a local four-quark operator between i?-meson states. 
There exists a vast literature on calculations of the decay consta.rit fp and the Bb 
parameter. CoT rihin ing the results of some recent QCD sum-rul(£3lir 
calculationsli23^E2^, we quote the value 



and lattice 



/s = 185±40 MeV. (186) 

Together with the prediction Bb — 1.08 obtained from lattice calculationf^H^, this 
gives 

B]^^fB ^ (200 ± 40) MeV . (187) 



Solving then (185) for \Vtd\, one obtains 



201 



\Vtd\ 



= 8.53 X 10- C-^] ^ITOGeVy - ( ^tub 



\ B^^^/b / V fntimt) J \0.465 ps ^ 



1/2 



(188) 



Taking mt{mt) = (170 ± 15) OeV for the running top-quark mass, and using the 
average experimental valualZ; 



Ams = (0.465 ±0.024) ps"^ , 



(189) 



gives 

|Vtd| = (8.53 ± 1.81) X 10^^ 
The corresponding range of values for Rt is 



Vtd 



0.99 ±0.22. 



(190) 



(191) 



6 






d 







Figure 18: Box diagrams for B^~B^ mixing in the Standard Model. 



Equations (184) and (191 



yield constraints on the Wolfenstein parameters p 
and fj, which have the form of rings centred at (p, ^) = (0,0) and (0,1). Another 
constraint can be obtained from the measurement of indirect CP violation in the 
kaon system. The experimental result on the parameter measuring CP violation 
in K^-K^ mixing implies that the unitarity triangle lies in the upper half plane. 
The constraint arising in the p-77 plane has the form of a hyperbola, the shape of 
which depends on a hadronic parameter Bk- The theoretical prediction isEy 



mt 
raw 



1.52 



+ (0.69 ±0.05) 



A^B 



K 



0.50, 



(192) 
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where A — 0.80 ± 0.04 according to (179). In the last few years, theoretical cal- 
culations of the Bk parameter have converged and give results in the ball park 
of 



Bk = 0.75 ±0.15. 



(193) 



In particular, the most recent lattice calculations 
the results obtained using the 1/Nc expansionr°^'P°^ 



12031 g^j-g gQQ(j agreement with 
and the differences yith pre- 
have 



vious, lower predictions for Bk based on duality and chiral symmetryE2£ 
been understood. 

In principle, the measurement of the ratio Re {e' /e) in the kaon system could 
provide a determinatio n of 77 i ndependent of p. In practice, howev e r , th e experimen- 



ever the 
2io|-pi2| 



tal situation is unclearE23Ei^!3, and the theoretical calculationsE£I^Ei3 of this ratio 
are affected by large uncertainties, so that there currently is no useful bound to be 
derived. 




Figure 19: Experimental constraints on the unitarity triangle in the p- 
fj plane. The region between the solid (dashed) circles is allowed by the 
measurement of Rb {Rt) discussed above. The dash-dotted curves show 
the constraint following from the measurement of the ek parameter in the 
kaon system. The shaded region shows the allowed range for the tip of the 
unitarity triangle. The base of the triangle has the coordinates (0, 0) and 
(1,0). 



In Fig. we show the constraints which the measurements of Rb, Rt, and ex 
imply in the p-fj plane. Given the present theoretical and experimental uncertainties 
in the analysis of charmless B decays, B^-B^ mixing, and CP violation in the kaon 
system, there is still a rather large region allowed for the Wolfenstein parameters. 
This has important implications. For instance, the allowed region for the angle (3 
of the unitarity triangle (see Fig. ^ is such that 



0.34 < sin 2/3 < 0.75 



(194) 
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Below, we will discuss that the CP asymmetry in the decay B i^Ks, which 
is one of the favoured modes to search for CP violation at a future B factory, is 
proportional to sin 2/3. Obviously, the prospects for discovering CP violation with 
such a machine depend on whether sin 2/3 is closer to the upper or lower bound in 



(194). A more reliable determination of the shape of the unitarity triangle is thus 
of great importance. 

On the other hand, our knowledge of the unitarity triangle has already improved 
a lot in the last few years, and we are now reaching a state where the analysis 
described in this section becomes a serious test of the Standard Model. If the three 
bands in Fig. |9| did not overlap, this would be an indication of New Physics. 

5.7. CP Asymmetries in Neutral B-Meson Decays 

As mentioned above, decays of neutral B mesons into CP eigenstates provide for 
largely model-independent determinations of the angles of the unitarity triangle. In 
the i?-meson system, up to corrections of order 10~^, we have 

This combination of CKM parameters can be read off directly from the vertices of 
the box diagrams in Fig. ^ which in the Standard Model are responsible for the 
non-diagonal element M*2 of the mass matrix. Notice that for the real part of the 
box diagrams, which determines M12, the contributions of c and u quarks in the 
loops can be neglected. 

To eliminate hadronic uncertainties, one has to choose decay modes dominated 
by a single diagram. However, most channels receive contributions from "tree" 



and 
ratioSi: 



41 



iiguin" diagrams, which for a generic b — > qq'q' decay contribute in the 
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penguin ^ ^tbVfij 

tree 12^ ' ' « 



ii^zri-'TT^- (196) 



The first factor arises from the loop suppression of the penguin diagrams and is 
of order 2%, the second factor accounts for the fact that hadronic matrix elements 
of penguin operators are usually enhanced with respect to those of the operators 
appearing in tree diagrams by r ~ 2-5, and the last factor is the ratio of CKM 
matrix elements. 

It follows that there are three possibilities to obtain the dominance of a single 
diagram: 

• If the CKM parameters of the penguin diagram are not enhanced with respect 
to the tree diagram, i.e. if 

VttV,; 



the tree diagram dominates over the penguin diagram. Examples of such 
decays are B tttt, B DD, Bg pKg and Bg i'Ks- 
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If tree diagrams are forbidden, the penguin diagram dominates. Examples of 
such decays are B <pKs, B KsKs, Bg — > rj'rj' and Bg 4>Ks- 



• If 



arg 



VtbVt 



tq 



= or TT , 



(198) 



both the tree and the penguin diagram have the same weak phase. In that 
case one stiU has |^/^| = 1, i.e. no hadronic uncertainties. Examples of such 
decays are B 4>Ks and Bs — > V"/'- 

Let us illustrate these three classes of decays with explicit examples. 



5.7.1. Tree-dominant decays: B — > tttt 

The decay B — > tttt proceeds through the quark decay b — > uud^ for which both 
the tree and the penguin diagram have CKM parameters of order A'^, as shown in 
Fig. BG. Thus, the tree diagram is dominant, and to a good approximation 



A -'^.^ 
p A 



V*uVud 



' td 



and 



Im A^TT — sin 2a . 



(199) 



(200) 



Hadronic uncertainties arise from the s ma ll admixture of penguin contributions, 
which lead to |A| ^ 1. They are exp ectedl^^^ to be of order 10%, and can be reduced 



further by using isospin analysis 



168 



172 





Figure 20: Tree and penguin diagrams for the decay B tttt. 



5.7.2. Tree-forbidden decays: B — > (f>Ks 

The decay B 4>Ks proceeds through the quark transition b sss, i.e. it involves a 
flavour-changing neutral current, which is forbidden at the tree level in the Standard 
Model. Thus, the relevant diagram is the penguin transition shown in Fig. |2l|. A 
new ingredient is the presence of K-K mixing, which adds a factor 



(201) 
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in the definition of A. This is essential for decays with a single Ks, since only 
and K'^ transitions are allowed, and interference between them is 

possible only due to K-K mixing. It follows that 



and therefore 



q 

PJ B 



pJk A 
ImA 



v^^ytd VcsV:^ VtbV; 



vt,v:, v:,vts 



sin 2^. 



(202) 



(203) 




Figure 21: Penguin diagram for the decay B 4>Ks. 



5.7.3. Decays with a single weak phase: B — > (fiKs 

The decay B ipKs is based on the quark transition h —t ccs, for which the tree 



diagram is dominant. As shown in Fig. 22, the tree amplitude is proportional to 



VrhV* - A2. One finds 



and therefore 



Im A^Ks ~ sin 2/3 . (205) 



In th e pr esent case, the contamination from the penguin contribution is extremely 



smalltZ^. Depending on the flavour q of the quark in the loop, the penguin con- 
tributions are proportional to VthVj* ~ A'^ (for q = t), VctV*^ ~ A^ (for q = c), 
and VubV*^ ~ A"^ (for q = u). Because of the relation VttV*^ = -VcbV*, + 0(A4), it 
follows that up to very small corrections the penguin contributions have the same 
weak phase as the tree diagram. Hence, their presence affects neither |A| nor ImA. 
Detailed estimates show that the hadronic uncertainties are only of order 10^'^. 
This makes the measurement of sin 2/3 in i? ^ "^Ks the theoretically cleanest de- 
termination of any CKM parameter. For this reason, this decay is often considered 
the "gold-plated" mode of a i? factory. 

The above-mentioned examples are only meant to illustrate the range of pos- 
sibilities for performing model-independent measurements of CP-violating CKM 
parameters in neutral _B-meson decays into CP eigenstates. A summary and some 
more examples are given in Table ^. The angle /3' appearing in the CP asymmetries 
for _Bs-meson decays is the analogue of the angle /3 in the unitarity triangle defined 
by the relation 

K.Ki + v^csK; + v^.i^;fc = o. (206) 
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be d d 




Figure 22: Tree and penguin diagrams for the decay B ip^S- 



Table 3: Examples of CP asymmetries for B and Bg decays into CP 
eigenstates 





B decays 


Bs decays 


Quark decay 


Final state 


SM prediction 


Final state 


SM prediction 


b — * CCS 




-sin 2/3 




-sin 2/3' 


b ccd 


D+D- 


-sin 2/3 




-sin 2/3' 


b — > uud 


TT+TT" 


sin 2a 


pKs 


-sin2(7 + /3') 


b sss 


<pKs 


-sin2(/3-/3') 


r]'r)' 





b ssd 


KsKs 





cf>Ks 


sin2(/3-/3') 



Experimentally, |sin2/3'| < 0.06. 

At the end of this section, let us stress again that the Standard Model descrip- 
tion of CP violation is at the same time very predictive (since all CP violation is 
related to a single parameter) and most likely wrong (because of the problems with 
baryogenesis and strong CP violation). Thus, the prospects arc good that once the 
various CP asymmetries in i?-meson decays can be explored at a i? factory, devi- 
ations from the picture described here will arise. Those deviations would indicate 
New Physics beyond the Standard Model. 

6. Concluding Remarks 

We have presented a review of the theory and phenomenology of heavy-flavour 
physics. The theoretical tools that allow to perform quantitative calculations in 
this area are the heavy-quark symmetry, the heavy-quark effective theory, and the 
l/rriQ expansion. We have discussed in detail exclusive weak decays of B mesons, 
inclusive decay rates and lifetimes of b hadrons, and CP violation in B-meson decays. 
Besides presenting the status of the latest developments in these fields, our hope 
was to convinced the reader that heavy-flavour physics is a rich and diverse area 
of research, which is at present characterized by a fruitful interplay between theory 
and experiments. This has led to many significant discoveries and developments on 
both sides. Heavy-quark physics has the potential to determine many important 
parameters of the electroweak theory and to test the Standard Model at low energies. 
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At the same time, it provides an ideal laboratory to study the nature of non- 
perturbative phenomena in QCD, still one of the least understood properties of 
the Standard Model. The phenomenon of CP violation, finally, is one of the most 
intriguing aspects of high-energy physics. Today, there is only a single unambiguous 
measurement of a CP-violating quantity. But already in a few years, when CP 
violation in the iJ-meson system can be explored at the B factories, this will very 
likely provide some clues about the physics beyond the Standard Model. 

Indeed, the prospects for further significant developments in the field of heavy- 
flavour physics look rather promising. With the approval of the first asymmetric 
B factories at SLAC and KEK, with ongoing i?- physics programs at the existing 
facilities at Cornell, Fermilab and CERN, and with plans for future B physics at 
HERA-B and the LHC-B, there are Beautiful times ahead of us! 
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